AFRL-OSR-VA-TR-201 5-0026 


ADVANCED  TECHNOLOGIES  FOR  STRUCTURAL  AND  FUCNTIONAL  OPITCAL  COHERENCE 
TOMOGRAPHY 


James  Fujimoto 

MASSACHUSETTS  INSTITUTE  OF  TECHNOLOGY 


01/07/2015 
Final  Report 


DISTRIBUTION  A:  Distribution  approved  for  public  release. 


Air  Force  Research  Laboratory 
AF  Office  Of  Scientific  Research  (AFOSR)/  RTB 
Arlington,  Virginia  22203 
Air  Force  Materiel  Command 


REPORT  DOCUMENTATION  PAGE 

Form  Approved 

OMB  No.  0704-0188 

Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and  maintaining  the 
data  needed,  and  completing  and  reviewing  this  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information,  including  suggestions  for  reducing 
this  burden  to  Department  of  Defense,  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports  (0704-0188),  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington,  VA  22202- 
4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  any  penalty  for  failing  to  comply  with  a  collection  of  information  if  it  does  not  display  a  currently 
valid  OMB  control  number.  PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS. 

1 .  REPORT  DATE  (DD-MM-YYYY)  2.  REPORT  TYPE 

1/2/2015  FINAL  REPORT 

3.  DATES  COVERED  (From  -  To) 

30  SEP  2010  to  29  SEP  2014 

4.  TITLE  AND  SUBTITLE 

ADVANCED  TECHNOLOGIES  FOR  STRUCTURAL  AND 

FUCNTIONAL  OPITCAL  COHERENCE  TOMOGRAPHY 

5a.  CONTRACT  NUMBER 

FA9550-1 0-1 -0551 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

6.  AUTHOR(S) 

JAMES  G.  FUJIMOTO 

Department  of  Electrical  Engineering  and  Computer  Science 

Research  Laboratory  of  Electronics 

Massachusetts  Institute  of  Technology 

50  Vassar  St.  Office  36-345 

5d.  PROJECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Massachusetts  Institute  of  Technology 

77  Massachusetts  Ave 

Cambridge,  MA  02139 

8.  PERFORMING  ORGANIZATION  REPORT 
NUMBER 

9.  SPONSORING  /  MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

HOWARD  SCHLOSSBERG  and  WILLIAM  P.  ROACH 

AFOSR/RTB 

875  North  Randolph  Road 

Suite  325,  Room  31 1 2 

Arlinnton.  VA  22203-1768 

10.  SPONSOR/MONITOR’S  ACRONYM(S) 

11.  SPONSOR/MONITOR’S  REPORT 
NUMBER(S) 

12.  DISTRIBUTION  /  AVAILABILITY  STATEMENT 

13.  SUPPLEMENTARY  NOTES 

14.  ABSTRACT 

This  program  is  a  collaborative  effort  between  investigators  at  the  Massachusetts  Institute  of  Technology;  the  New  England  Eye 
Center,  Tufts  University  School  of  Medicine;  and  the  Martinos  Center  for  Biomedical  Imaging,  Massachusetts  General  Hospital. 
The  program  focuses  on  the  development  and  application  of  advanced  technology  for  OCT  with  an  emphasis  on  ultrahigh 
speed  structural  and  functional  imaging.  Emphasis  is  placed  on  applications  to  ophthalmology,  the  area  in  which  OCT  has  had 
the  largest  clinical  impact,  and  brain  science,  an  emerging  area  where  structural  and  functional  imaging  are  enabling  advances 
in  fundamental  research.  Functional  imaging  of  the  retina  and  brain  are  synergistic  and  provide  a  cross-platform  approach  to 
develop  and  validate  new  methods.  The  advanced  technology  and  methods  developed  in  this  program  promise  to  enable  new 
applications  in  multiple  areas  of  biomedical  research  and  clinical  medicine. 

15.  SUBJECT  TERMS 

Biomedical  optics,  optical  imaging,  optical  coherence  tomography,  ophthalmology,  neurosciences,  brain  imaging 

16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION 

OF  ABSTRACT 

18.  NUMBER 
OF  PAGES 

19a.  NAME  OF  RESPONSIBLE  PERSON 

a.  REPORT 

b.  ABSTRACT 

c.  THIS  PAGE 

19b.  TELEPHONE  NUMBER  (include  area 
code) 

Standard  Form  298  (Rev.  8-98) 


Prescribed  by  ANSI  Std.  Z39.18 


AFOSR  Grant  FA9550-10-1-0551 

Advanced  Technologies  for  Structural  and  Functional 
Optical  Coherence  Tomography 
September  30,  2010  to  September  29,  2014 
Final  Report 


Table  of  Contents 

1 .0  Introduction  2 

1 . 1  Program  Highlights  3 

2.0  Next  Generation  OCT  Technology  4 

2.1  MEMS-tunable  VCSELs  for  Ultrahigh  Speed  OCT  4 

2.2  Functional  OCT  Imaging  5 

OCT  Angiography  (OCTA)  5 

Doppler  OCT  6 

Polarization  Sensitive  OCT  (PS-OCT)  7 

2.3  Ultralong  Range  Imaging  using  MEMS  VCSEL  Technology  10 

2.4  Development  of  Motion  Correction  Algorithms  for  Ophthalmology  1 1 

3.0  Ophthalmic  OCT  13 

3.1  Technology  Development  for  Clinical  Ophthalmology  14 

3.2  Handheld  OCT  16 

3.3  Wide  Field  Retinal  Imaging  18 

3.4  Long  Range  Imaging  of  the  Anterior  Eye  and  Full  Eye  Length  Imaging  19 

3.5  OCT  Angiography  of  the  Retina  and  Choroid  21 

Current  Techniques  for  Visualizing  Retinal  and  Choriocapillaris  Microvasculature  2 1 
Development  of  SS-OCT  Angiography  using  the  VCSEL  Swept  Light  Source  2 1 

Demonstration  of  SS-OCT  Angiography  Imaging  in  Normal  Subjects  22 

Exudative  AMD  23 

Dry  AMD  23 

Diabetic  Retinopathy  25 

3.6  En  face  Doppler  OCT  Measurement  of  Total  Retinal  Blood  Flow  26 

Measurement  of  Total  Retinal  Blood  Flow  26 

Conventional  Techniques  for  TRBF  Measurement  26 

En  Face  Doppler  OCT  27 

En  Face  Doppler  OCT  with  High  Speed  OCT  Technology  27 

Synchronized  En  Face  Doppler  OCT  with  Pulse  Oximetry  Cardiac  Gating  28 

Clinical  Studies  of  Total  Retinal  Blood  Flow  28 

4.0  Brain  Imaging  29 

4. 1  Introduction  to  OCT  for  Neuroscience  Applications  29 

4.2  Quantitative  Measurements  of  Absolute  Cerebral  Blood  Flow  (CBF)  3 1 

4.3  DLS-OCT  Imaging  of  CBF  and  Cellular  Viability  3 1 

4.4  OCT  Imaging  of  Capillary  RBC  Flow  Properties  33 

4.4  Rapid  Volumetric  Imaging  of  Capillary  Network  Flux  34 


1 


5.0  Endoscopic  OCT  Technology  and  Applications  36 

5.1  Clinical  Studies  with  Endoscopic  OCT  36 

Characterization  of  Buried  Glands  Pre-  and  Post-Radiofrequency  Ablation  with  OCT  37 
Endoscopic  OCT  Markers  of  Radiofrequency  Ablation  Treatment  Response  39 

5.2  Endoscopic  OCT  Technology  Development  40 

Ultrahigh  Speed  Endoscopic  Swept  Source  OCT  40 

Micromotor  Endoscopic  Imaging  Probes  41 

5.3  Novel  Diagnostics  with  Endoscopic  OCT  42 

En  Face  OCT  and  Pitt  Patterns  43 

Endoscopic  OCT  A  and  Microvasculature  45 

6.0  Pathology  Imaging  47 

6. 1  Multiphoton  Microscopy  for  Rapid  Evaluation  of  Pathology  48 

Introduction  to  multiphoton  microscopy  (MPM)  48 

Preliminary  Multiphoton  Microscopy  Study  of  Breast  Surgical  Specimens  49 

Development  of  Ultrahigh  Speed  Wide  Field  of  View  Multiphoton  Microscopy 

for  Pathology  5 1 

Development  of  Gigavoxel  Multiphoton  Fluorescent  Lifetime  Imaging  Technology  52 

6.2  Development  of  Next  Generation  Optical  Coherence  Microscopy  Technology  54 

Optical  Coherence  Microscopy  (OCM)  54 

High  Speed  Fourier  Domain  Optical  Coherence  Microscopy  55 

Image  Tiling  and  Display  for  Large  Field  of  View  57 

7.0  Other  Collaborations  58 

7.1  Computer  Aided  Diagnosis  of  Plaque  Erosion  58 

7.2  OCT  for  Evaluation  of  Kidney  Transplantation  61 

8.0  Publications  under  AFOSR  program  62 


1.0  Introduction 

This  program,  AFOSR  FA9550-10-1-0551,  entitled  “Advanced  Technologies  for  Structural  and 
Functional  Optical  Coherence  Tomography,”  is  a  multidisciplinary  collaborative  program  among 
investigators  at  the  Massachusetts  Institute  of  Technology,  Department  of  Electrical  Engineering 
and  Computer  Science  and  Research  Laboratory  of  Electronics  (P.I.  James  Fujimoto,  Elihu 
Thomson  Professor  of  Electrical  Engineering  and  Computer  Science);  the  Harvard  Medical 
School,  Massachusetts  General  Hospital,  Martinos  Center  for  Biomedical  Imaging  (David  Boas, 
Professor  in  Radiology  at  Harvard  Medical  School,  Director,  Optical  Imaging  Core  &  Lab, 
Martinos  Center,  Department  of  Radiology,  MGH)  and  the  Tufts  University  School  of  Medicine, 
New  England  Eye  Center  (Jay  Duker,  M.D.,  Director  of  the  New  England  Eye  Center,  Chairman 
of  Ophthalmology).  The  objective  of  this  program  is  to  develop  the  next  generation  of  advanced 
technology  for  OCT  with  an  emphasis  on  ultrahigh  speed  structural  and  functional  imaging  and 
to  explore  applications  in  key  areas.  Emphasis  is  placed  on  ophthalmology,  the  area  in  which 
OCT  has  had  the  largest  clinical  impact  to  date,  and  brain  science,  an  emerging  area  where 
structural  and  functional  imaging  are  enabling  advances  in  fundamental  research.  Functional 
imaging  of  the  retina  and  brain  are  synergistic  and  provide  a  cross-platform  approach  to  develop 
and  validate  new  methods.  The  technology  and  approaches  developed  in  this  program  promise  to 
lay  the  foundation  for  a  range  of  new  applications  in  fundamental  research  and  clinical  medicine. 
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1.1  Program  Highlights 


This  program  was  in  a  one  year  no  cost  extension  which  ended  in  September  2014  and  therefore 
this  Final  Report  describes  results  over  a  four  year  period.  The  highlights  of  the  program  include 
technological,  scientific  and  clinical  advances,  as  well  as  technology  transfer  and  new  program 
development.  New  technology  was  developed  and  translated  for  studies  in  multiple  clinical 
applications  in  ophthalmology  and  endoscopy.  Technology  and  intellectual  property  developed 
under  the  program  was  transferred  to  industry  and  released  a  commercial  ophthalmic  imaging 
product.  AFOSR  funding  also  supported  pilot  studies  in  new  areas  which  subsequently  lead  to 
the  development  of  new  NIH  funded  research  programs.  Highlights  of  major  accomplishments  in 
this  program  include: 

•  Publication  of  over  75  peer  reviewed  journal  articles. 

•  Collaboration  with  Thorlabs  /  Praevium  Research  to  develop  MEMS-VCSEL  lasers  as  a 
platform  technology  for  structural  and  functional  swept  source  OCT  imaging. 
Demonstration  of  multiple  new  ophthalmic  and  endoscopic  imaging  applications. 

•  Invention  of  software  based  motion  correction  with  collaborators  at  the  University  of 
Erlangen.  The  technology  was  transferred  to  industry  and  released  in  2014  an  ophthalmic 
imaging  instrument  (Avanti  by  Optovue)  outside  of  the  US  and  is  pending  clearance  in 
the  US. 

•  Development  and  demonstration  of  handheld  ophthalmic  OCT  technology  for  3D  retinal 
imaging. 

•  Development  of  prototype  OCT  technology  for  ophthalmology  which  is  used  for  clinical 
studies  at  the  New  England  Eye  Center,  the  University  of  Pittsburgh  Eye  Center  and  the 
Casey  Eye  Institute,  Oregon  Health  and  Sciences  University. 

•  Development  of  new  high  performance  OCT  technology  for  endoscopic  imaging  which  is 
used  for  clinical  studies  at  the  Boston  VA  Healthcare  System,  Harvard  Medical  School. 

•  OCT  angiography  studies  of  the  choriocapillaris  in  exudative  AMD.  This  study  was 
featured  on  the  cover  of  Ophthalmic  Surgery  Lasers  and  Imaging  Retina,  45,  6,  Nov/Dec. 
2014,  Moult,  et  al.  [1] 

•  First  demonstration  of  endoscopic  OCT  angiography  enabling  three  dimensional  imaging 
of  subsurface  vasculature  in  the  GI  tract.  This  study  was  featured  on  cover  of 
Gastroenterology,  147,  6,  December  2014,  Tsai,  et  al.  [2] 

•  Demonstration  of  nonlinear  microscopy  for  assessing  cancer  pathology  in  breast  surgical 
specimens.  Published  in  PNAS,  1 1 1,  43,  2014,  Tao,  et  al.  [3],  this  study  demonstrated 
that  nonlinear  microscopy  can  assess  cancer  pathology  in  surgical  specimens  with 
comparable  sensitivity  and  specificity  to  standard  histology  and  promises  to  reduce  rates 
of  second  surgeries. 

•  Development  of  a  new  NIH  program  (R01  CA178636-02,  “Intraoperative  real  time  breast 
cancer  margin  assessment  with  nonlinear  microscopy”,  P.I.  James  Fujimoto)  in 
collaboration  with  investigators  at  the  Harvard  Medical  School  and  Beth  Israel 
Deaconess  Medical  Center. 

•  Development  of  a  new  NIH  program  (R01  DK094877  -01,  “Non-invasive  evaluation  of 
transplant  kidney  using  OCT”,  P.I.  Peter  Andrews)  in  collaboration  with  investigators  at 
the  Georgetown  University  and  University  of  Maryland. 
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2.0  Next  Generation  OCT  Technology 

2.1  MEMS-tunable  VCSELs  for  Ultrahigh  Speed  OCT 

Optical  coherence  tomography  (OCT)  is  a  noninvasive  imaging  modality  which  can  visualize 
cross-sectional  and  three-dimensional  tissue  microstructure  with  high  resolution  and  imaging 
speeds.  OCT  is  analogous  to  ultrasound,  measuring  the  depth-dependent  backscattered  power  of 
infrared  light  instead  of  acoustic  waves.  OCT  is  a  powerful  technology  for  biomedical  research 
and  clinical  diagnostics  because  tissue  can  be  examined  noninvasively  in  real-time,  without 
excisional  biopsy  and  processing  required  for  histopathology.  OCT  has  excellent  imaging 
resolution  (2  to  10  microns)  which  enables  visualization  of  tissue  structure  and  pathology.  OCT 
has  been  demonstrated  using  multiple  patient  interfaces  and  medical  devices  including 
ophthalmoscopes,  endoscopes,  catheters,  needles,  and  other  clinical  or  surgical  instruments. 
Extensions  of  OCT  technology  including  Doppler  OCT,  OCT  angiography  (OCTA)  and 
polarization  sensitive  OCT  (PS-OCT)  enable  integrated  structural  and  functional  imaging. 

A  major  focus  of  this  program  has  been  on  developing  advanced  OCT  technology  for 
structural  and  functional  imaging,  with  the  aim  of  creating  new  methods  which  would  be 
applicable  across  a  wide  range  of  biomedical  applications  and  clinical  specialties.  Recent 
advances  in  wavelength-swept  lasers  enable  ultrahigh-speed  imaging,  opening  up  a  wide  variety 
of  functional  OCT  imaging  applications.  During  the  previous  contract  period,  our  group  at  MIT 
worked  in  collaboration  with  Praevium  Research  and  Thorlabs  to  develop  MEMS-tunable 
vertical  cavity  surface  emitting  lasers  (VCSELs)  operating  at  1310  mn  and  1050  mn  center 
wavelengths  and  demonstrate  swept  source/Fourier  domain  OCT  (SS-OCT)  in  key  ophthalmic 
and  endoscopic  applications.  Figure  1  shows  a  schematic  diagram  and  photograph  of  the  MEMS 
VCSEL  laser.  The  MEMS-tunable  VCSELs  are  wavelength- swept,  semiconductor  lasers  that 
emit  a  beam  perpendicular  to  the  mirror  surface.  One  of  the  key  features  of  the  MEMS-tunable 
VCSEL  technology  is  the  micro-scale  cavity  length  and  high  MEMS  resonance  frequency  that 
allow  single  longitudinal  mode,  ultra-narrow  linewidth  operation  at  ultrahigh  laser  sweep  rates 
(>1  MHz).  This  corresponds  to  15  to  40-fold  increase  in  imaging  speed  compared  to 
commercially  available  spectral  /  Fourier  domain  (SD-OCT)  systems.  Using  this  new  technology, 
our  collaborative  group  demonstrated  ophthalmic  imaging  at  400  kHz  to  580  kHz  A-scan  rates  [4, 
5]  as  well  as  endoscopic  imaging  at  1  MHz  A-scan  rate  [6]. 

Swept  source  OCT  using  the  VCSEL  laser  enables  both  high  speed  as  well  as  long  imaging 
range  and  has  been  a  key  technology  in  this  program.  High  speed  is  especially  powerful  because 
it  enables  functional  imaging  such  as  Doppler  flow  measurement,  motion  contrast  using  OCT 
angiography  (OCTA),  or  polarization  sensitivity  OCT  (PS-OCT). 
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Figure  1.  Schematic  and  photograph  of  a  VCSEL  laser.  The  laser  consists  of  a  semiconductor  mirror,  quantum 
well  gain  material  and  dielectric  mirror  fabricated  on  a  MEMS.  The  gain  is  optically  pumped  by  a  laser  diode 
(not  shown).  The  VCSEL  cavity  is  short  and  can  operate  in  a  single  longitudinal  mode,  providing  very  narrow 
linewidth.  The  MEMS  mirror  is  electrostatically  actuated  and  can  vary  the  cavity  length  to  sweep  the  laser 
frequency  /  wavelength.  The  MEMS  mirror  can  achieve  rapid  sweep  rates  as  well  as  adjustable  sweep  ranges. 


2.2  Functional  OCT  Imaging 

OCT  angiography  (OCTA) 

OCT  angiography  (OCTA)  is  a  functional  OCT  imaging  technique  that  enables  three 
dimensional  visualization  of  microvasculature  without  requiring  intravenous  dyes,  such  as 
fluorescein  or  indocyanine  green.  OCTA  is  especially  powerful  for  ophthalmology  research  and 
clinical  investigation.  Retinal  and  choroidal  microvasculature  are  imaged  by  measuring 
fluctuations  caused  by  the  passage  of  red  blood  cells  through  an  OCT  voxel  by  rapidly  acquiring 
repeated  B-scans  from  the  same  location  [7-16].  Because  the  same  location  must  be  re-scanned 
multiple  times,  volumetric  OCTA  imaging  requires  much  larger  numbers  of  A-scans  than 
structural  OCT  and  ultrahigh  imaging  speeds  are  important.  Although  OCT  angiography  cannot 
visualize  vascular  leakage  as  is  possible  with  conventional  angiography,  OCTA  can  be 
performed  repeatedly  because  exogenous  contrast  agents  are  not  administered.  Furthermore, 
because  OCT  angiography  is  based  on  motion  contrast,  it  can  potentially  provide  information  on 
local  retinal  blood  perfusion  at  the  capillary  level  in  addition  to  three  dimensional  structural 
information. 

Figures  2  and  3  show  retinal  and  choroidal  OCT  angiograms  acquired  using  a  1050  nm 
wavelength,  ultrahigh  speed  swept  source  OCT  prototype  instrument  operating  at  a  400  kHz  A- 
scan  rate.  Figure  2  shows  OCT  retinal  angiograms  of  12  mm  x  12  mm  and  3  mm  x  3  mm 
regions  acquired  from  a  normal  volunteer.  Although  the  large  field  of  view  OCT  angiogram  in 
Figure  2(A)  visualizes  large  vascular  structure,  visualization  of  progressively  finer  details 
requires  scanning  smaller  fields  of  view  as  in  Figure  2(B).  Figure  3  shows  a  wide-field 
mosaicked  OCTA  spanning  -32  mm  on  the  retina  at  the  level  of  choriocapillaris,  the  capillary 
layer  of  the  choroid,  acquired  from  a  normal  volunteer.  Multiple  3  mm  x  3  mm  field  of  view 
OCTAs  were  stitched  together  to  generate  a  panoramic  visualization  of  the  choriocapillaris, 
demonstrating  variations  in  the  choriocapillaris  density  and  structure  at  different  locations, 
consistent  with  histological  studies  [17].  Ultrahigh  speed  is  crucial  for  OCT  angiography  because 
it  enables  larger  retinal  coverage  without  stitching  and/or  higher  resolution  visualization  of 
microvasculature.  Studies  of  OCTA  in  ophthalmic  and  endoscopic  applications  are  further 
described  in  sections  3.5  and  5.3. 
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Figure  2.  OCT  angiograms  (OCTA)  of  retinal  capillary  structure.  (A)  12  mm  x  12  mm  and  (B)  3  mm  x  3 
mm  fields  showing  wide  field  retinal  vasculature  as  well  as  the  foveal  avascular  zone.  Capillary  structure 
can  be  imaged  in  3D  and  image  analysis  can  be  used  to  quantitatively  analyze  the  vascular  network. 
Exogenous  contrast  is  not  required.  Scale  bars:  1  mm. 


Figure  3.  OCT  angiograms  (OCTA)  of  choriocapillaris  microvasculature  spanning  ~32  mm  on  the  retina. 

(A)  Mosaicked  OCT  intensity  en  face  projection  images.  (B)  Mosaicked  OCT  angiograms  of  the 
choriocapillaris.  Different  microvasculature  patterns  and  densities  can  be  observed  at  different  fundus 
locations,  consistent  with  known  histological  findings.  Scale  bars:  1.5  mm. 

Doppler  OCT 

Doppler  OCT  is  another  functional  OCT  imaging  technique  which  is  analogous  to  ultrasound 
color  Doppler  imaging.  Doppler  OCT  detects  the  velocities  of  moving  scatterers  such  as 
erythrocytes  flowing  inside  blood  vessels.  Since  Fourier  domain  OCT  acquires  phase 
information  along  with  back-reflected  intensity,  Doppler  measurements  can  be  performed  by 
calculating  the  difference  in  optical  phase  between  subsequent  A-scans.  In  ophthalmic 
applications,  retinal  blood  flow  can  be  assessed  by  scanning  large  blood  vessels  in  or  around  the 
optic  disc.  High  imaging  speed  is  beneficial  for  Doppler  OCT  since  Doppler  requires  the  spatial 
A-scan  sampling  interval  to  be  much  closer  than  the  optical  spot  size,  so  that  A-scans  have 
correlated  speckle  information  and  phase  shifts  from  blood  flow  can  be  measured.  High  speed 
also  enables  new  measurement  techniques  such  as  en  face  Doppler,  where  blood  flow  is 
measured  in  en  face  planes  using  a  volumetric  data  set.  As  described  in  section  4,  this  technique 
was  originally  developed  by  Dr.  David  Boas’  group  at  MGH  for  measuring  blood  flow  in  the 
brain  and  was  translated  to  ophthalmology  [18].  En  face  Doppler  has  the  advantage  that 
determination  of  the  Doppler  angle  of  the  blood  vessel  is  not  required,  because  the  total  flow  can 
be  measured  by  integrating  the  velocity  over  the  vessel  area  in  the  en  face  plane.  However,  en 
face  Doppler  requires  extremely  high  imaging  speeds  because  multiple  volumes  must  be 
acquired  rapidly  enough  to  resolve  pulsatile  flow  from  the  cardiac  cycle.  Figure  4  shows  the  en 
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face  mapping  of  blood  flow  velocity  in  the  central  retinal  artery  at  different  phases  in  the  cardiac 
cycle.  Applications  of  Doppler  OCT  to  ophthalmology  iare  discussed  in  section  3.6. 


Figure  4.  En  face  Doppler  OCT.  (A)  En  face  OCT  of  human  optic  disc  and  (B,  C)  blood  flow  velocity  map 
measured  by  Doppler  OCT.  A  1.5  mm  x  2  mm  area  was  raster  scanned  to  generate  volumetric  OCT  data.  A 
blood  flow  velocity  map  was  extracted  from  the  automatically  selected  en  face  plane  at  an  optimal  depth. 

Units  pL/min,  scale  bars:  500  um. 

Polarization  sensitive  OCT  (PS-OCT) 

PS-OCT  is  a  functional  extension  of  conventional  OCT  that  can  assess  depth-resolved  tissue 
birefringence  in  addition  to  intensity  [19].  Examples  of  birefringent  tissue  include:  corneal 
stroma,  sclera  and  retina  nerve  fiber  layer  (RNFL),  collagen,  muscle,  tendons,  nerve,  bones, 
cartilage  and  teeth.  Alteration  of  tissue  birefringence  is  often  a  hallmark  for  disease  progression. 
PS-OCT  has  been  applied  for  anterior  and  posterior  eye  imaging  [20],  burn  depth  and  thermal 
damage  assessment  [21],  and  collagen  assessment  for  atherosclerotic  plaque  characterization 
[22].  Although  many  different  PS-OCT  techniques  have  been  demonstrated  [23-29],  their 
clinical  translation  has  been  complicated  by  their  high  degree  of  complexity.  Free  space  PS-OCT 
systems  are  susceptible  to  misalignment  and  are  difficult  to  operate  [25].  Fiber  based  PS-OCT 
using  polarization  maintaining  (PM)  fibers  requires  compensation  of  PM  fiber  mismatch 
between  reference  and  sample  arms  in  post-processing  [24].  Moreover,  these  PS-OCT  techniques 
require  known  polarization  states  of  the  incident  light  on  the  sample,  which  is  difficult  to  achieve 
in  catheter  or  endoscope-based  OCT.  Fiber-based  PS-OCT  using  polarization  multiplexing 
illuminates  the  sample  with  two  or  more  polarization  states  [27-29]  and  can  be  used  in  catheter- 
based  OCT.  Ffowever,  most  existing  fiber-based  PS-OCT  systems  are  complex  and  require  high 
system  SNR,  long  imaging  range  and  high  phase  stability. 

We  have  developed  a  simpler,  more  compact  and  robust  all-fiber  swept  source  PS-OCT 
technique  that  is  particularly  suitable  for  clinical  translation.  The  system  diagram  is  shown  in 
Figure  4.  Polarization  multiplexing  and  polarization  sensitive  detection  are  used  to  extract  tissue 
birefringence.  Polarization  multiplexing  is  performed  using  a  polarization  maintaining  (PM) 
fiber  in  the  sample  path.  PM  fiber  has  polarization  mode  dispersion  (PMD),  and  by  selecting  a 
proper  length  of  PM  fiber,  signals  from  the  two  orthogonal  polarization  states  of  the  input  light 
can  be  acquired  in  a  single  A-scan,  but  separated  at  different  depths.  For  polarization  sensitive 
detection,  two  fiber-based  polarization  beam  splitters  (FPBS)  are  used  to  eliminate  free  space  to 
fiber  coupling  and  improve  sensitivity.  Polarization  controllers  are  used  to  balance  the  FPBS 
based  polarization  sensitive  detection,  which  can  be  shown  by  Jones  Matrix  analysis  to 
unambiguously  resolve  the  tissue  birefringence. 
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Fiber-based  OCT  systems  greatly  augment  the  ability  of  OCT  by  enabling  fiber  coupled, 
handheld  OCT  instruments  as  well  as  catheter,  endoscope  and  needle  based  imaging  which  can 
access  luminal  organ  systems  such  as  the  gastrointestinal  and  pulmonary  tracts  as  well  as  solid 
organs  and  tumors  [30].  We  have  recently  demonstrated  a  catheter  /  endoscope  based  PS-OCT 
system  with  an  all-fiber  PS-OCT  design  (Figure  5)  using  a  MEMS-tunable  VCSEL  light  source 
[31]  and  a  micromotor  imaging  probe  [6].  Figure  6  shows  preliminary  imaging  results  from  the 
PS-OCT  system  in  human  finger,  lip,  ex  vivo  swine  esophagus  and  coronary  arteries.  It  can  be 
seen  that  PS-OCT  provides  tissue  birefringence  contrast  which  can  differentiate  layers  inside  the 
tissue.  For  example,  in  the  esophagus,  epithelium  is  associated  with  low  birefringence,  whereas 
lamia  propria  and  muscularis  mucosa  have  high  birefringence  (Figure  6  (F-G)). 


PM  fiber 


Figure  5.  Schematic  of  the  depth-encoded  all-fiber  PS-OCT  system.  The  system  uses  a  short  length  of  PM 
fiber  to  generate  two  time  multiplexed  incident  polarization  states.  The  OCT  interference  signal  is  detected 
with  a  polarization  sensitive  dual  balanced  detector.  This  technique  is  simpler  and  more  robust  than 
previous  approaches  and  promises  to  enable  a  wide  range  of  clinical  applications. 


Epi  dermis 


Intima 

Media 

Adventitia 


Figure  6.  Catheter-based  PS-OCT  imaging.  (A)  A  single  frame  OCT  intensity  image  of  thick  skin  from  a 
human  index  finger  tip.  (E)  Corresponding  phase  retardation  image  shows  high  retardance  at  the  dermal 
papillae.  (B)  A  single  frame  intensity  image  of  human  lip.  (V)  Corresponding  phase  retardation  image 
shows  low  retardance  at  the  epithelium  (EP)  and  high  retardance  at  the  lamia  propria  (LP).  (C)  Single 
frame  intensity  image  shows  the  epithelium  (EP)  and  lamia  propria  (LP)  /  muscularis  mucosa  (MM)  layers 
of  swine  esophagus  ex  vivo.  (G)  Corresponding  phase  retardation  image  shows  low  retardance  at  EP  but 
high  retardance  at  LP/MM.  (D)  Single  frame  intensity  image  of  a  swine  right  coronary  artery  ex  vivo.  (H) 
Corresponding  phase  retardation  image  shows  low  retardance  in  the  intima  layer,  and  high  retardance  in  the 
media  and  adventitia.  Horizontal  scale  bar:  30  degrees;  vertical  scale  bar:  500  um. 
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OCT  speckle  noise  can  significantly  affect  polarimetry  measurement  and  must  be  reduced  for 
birefringence  calculation  [32-34].  Conventional  speckle  reduction  methods  which  use  complex 
data  averaging  at  surface  Jones  matrices  [32],  or  directly  on  the  stokes  vectors  [33]  or  Jones 
matrices  [35],  require  phase  stability  over  successive  A-scans.  However,  for  SS-OCT,  laser 
sweep  jitter  or  sampling  fluctuation  can  make  the  phase  unstable.  Existing  phase  stabilization 
methods  complicate  PS-OCT  designs  [5,  36].  In  addition,  if  phase  stability  is  required,  it  is  not 
possible  to  average  data  over  multiple  acquisitions  or  data  which  has  been  motion  corrected 
using  automated  software  [37].  We  developed  a  simplified  speckle  reduction  method  which  does 
not  require  phase  stability,  thereby  greatly  simplifying  system  design  and  implementation.  By 
multiplying  the  measured  Jones  matrices  by  the  inverse  of  the  surface  Jones  matrix  [26],  we  can 
generate  a  matrix  similar  to  a  diagonal  matrix.  The  phase  difference  of  the  two  eigenvalues  of  the 
similar  matrix  contains  the  phase  retardation  information  from  the  tissue.  It  can  be  shown  that  the 
two  eigenvalues  suffer  from  the  same  unstable  phase.  By  sorting  the  eigenvalues  with  respect  to 
their  magnitude,  averaging  can  be  directly  performed  over  the  eigenvalue  ratio,  cancelling  the 
unstable  phase.  Preliminary  retinal  imaging  results  utilizing  averaging  over  a  small  kernel  and 
over  repeated  frames  are  shown  in  Figure  7.  This  technique  enables  more  accurate  polarimetry 
measurement  and  quantitative  assessment  of  tissue  birefringence. 


Figure  7.  Retinal  PS-OCT  imaging  with  speckle  noise  reduction.  The  scan  area  is  3  mm  x  3  mm  with  300 
x  300  x  5  A-lines/volume  (5  repeated  B-scans  per  location).  (A)  OCT  fundus  image.  (B)  Single  frame 
cross-sectional  intensity  image.  (C)  Corresponding  phase  retardation  image  without  speckle  noise  reduction 
(D)  Phase  retardation  image  using  speckle  noise  reduction.  Scale  bar:  500  um. 


Esophageal  cancer  is  one  of  the  most  lethal  malignancies,  with  a  five-year  survival  rate  of 
only  16%  [38].  Barrett’s  Esophagus  (BE)  is  associated  with  gastroesophageal  reflux  disease 
(GERD)  and  BE  with  dysplasia  confers  an  increased  risk  of  esophageal  adenocarcinoma  [39]. 
Endoscopic  RFA  is  a  therapy  [40-44]  utilizing  electrode  arrays  to  deliver  radiofrequency  (RF) 
energy  superficially  to  ablate  BE  with  a  low  stricture  rate  [45,  46].  However,  repeated  RFA 
treatments  are  required  to  achieve  complete  eradication  of  intestinal  metaplasia  (CE-IM)  [47-49]. 
Improving  the  effectiveness  of  RFA  would  reduce  the  number  of  treatments,  improve  cost 
effectiveness  and  make  therapy  more  widely  available.  However,  currently  there  are  no  effective 
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methods  to  directly  evaluate  the  efficacy  of  RFA.  Our  previous  studies  showed  that  endoscopic 
OCT  could  identify  markers  of  treatment  response  based  on  BE  thickness  pre-treatment  or 
residual  glandular  structure  post-treatment  [50].  In  addition  to  structural  OCT,  our  preliminary 
studies  also  suggest  that  RFA  results  in  loss  of  tissue  birefringence  due  to  thermal  damage  of  the 
lamia  propria,  which  can  be  assessed  by  PS-OCT.  Figure  8  shows  endoscopic  PS-OCT  imaging 
of  the  swine  esophagus  ex  vivo  before  Figure  8  (A,  B))  and  after  RFA  (Figure  8  (C,  D)).  Before 
RFA,  the  low  retardance  epithelium  can  be  clearly  differentiated  from  the  underlying  high 
retardance  lamina  propria  /  muscularis  mucosa  layer.  After  RFA,  the  low  retardance  region 
extended  deeper,  indicating  loss  of  birefringence  caused  by  the  ablation. 


Figure  8.  Endoscopic  PS-OCT  imaging  for  monitoring  RFA.  Preliminary  results  from  the  swine  esophagus 
ex  vivo.  RFA  was  applied  2x  with  an  energy  of  15  J/cm2,  consistent  with  clinical  treatment  protocols. 
Intensity  and  phase  retardation  images  of  swine  esophagus  before  (A,  B)  and  after  RFA  (C,  D).  The  phase 
retardation  images  were  averaged  over  10  frames  after  non-rigid  registration  to  correct  catheter  non- 
uniform  rotational  distortion.  Horizontal  scale  bar:  30  degrees;  vertical  scale  bar:  500  um. 


2.2  Ultralong  Range  Imaging  Using  MEMS  VCSEL  Technology 


The  micron-scale  cavity  length  of  MEMS-tunable  VCSELs  enables  single  longitudinal  mode 
operation,  resulting  in  an  extremely  long  coherence  length,  estimated  to  be  in  excess  of  100 
meters.  This  is  significantly  longer  than  the  coherence  length  of  other  swept  lasers  using  a  short 
cavity,  which  is  on  the  order  of  millimeters  to  1  cm.  Whereas  the  imaging  range  of  traditional 
SS-OCT  systems  has  been  limited  by  laser  coherence  length,  VCSEL  systems  are  only  limited 
by  the  analog  detection  bandwidth  and  A/D  sampling  rate,  potentially  enabling  micron-scale 
measurements  of  macroscopic  objects.  The  VCSEL  can  be  adjustably  scanned  with  a 
programmable  drive  waveform,  enabling  the  available  analog  bandwidth  to  be  optimized  for 
imaging  range,  axial  resolution,  or  A-scan  rate.  We  utilized  the  long  imaging  range  to  visualize 
the  anterior  segment  of  the  human  eye  which  comprises  the  cornea,  iris  and  crystalline  lens  [4], 
The  results  are  summarized  in  section  3.4. 

In  addition  to  biometry  and  full-eye  imaging,  the  ability  to  image  many  centimeter  scale 
objects  at  micron  resolution  may  have  important  industrial  or  military  applications  in 
nondestructive  testing  and  validation.  To  demonstrate  this  capability,  we  performed  long  range 
OCT  imaging  of  large  objects,  including  a  high  aspect  ratio  machined  metal  optical  housing, 
light  bulbs,  and  ammunition  [51].  Figure  9  shows  examples  of  long  range  OCT  imaging. 
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Figure  9.  Applications  of  long  range  OCT  imaging.  Imaging  examples  of  (A)  Optical  post  holder  which 
has  a  high  aspect  ratio  bore  hole.  (B)  Bullets  in  colored  gelatin.  (C)  Light  bulb.  White  arrow  indicates  the 
filament  holder. 


2.4  Development  of  Motion  Correction  Algorithms  for  Ophthalmology 

Working  in  collaboration  with  Prof.  Joachim  Hornegger’s  group  at  the  Friedrich  Alexander 
University  of  Erlangen-Nuremberg  in  Germany,  we  have  developed  new  software  based  motion 
correction  technology  for  ophthalmic  OCT  that  corrects  eye  motion  in  transverse  and  axial 
directions  and  enables  multiple  volumetric  data  sets  to  be  merged  to  increase  signal  to  noise  [37, 
52],  Intellectual  property  was  fded  based  on  this  research  (US  Patent  Application  13/097785, 
Filed  April  29,  2011,  Kraus,  et  al.  “Method  and  Apparatus  for  Motion  Correction  and  Image 
Enhancement  for  Optical  Coherence  Tomography”).  The  technology  was  licensed  and 
transferred  to  Optovue,  Inc.  and  released  in  a  commercial  ophthalmic  OCT  instrument  (the 
Optovue,  RTVue-XR)  in  2014,  which  is  available  outside  the  United  States.  Approval  in  the  US 
is  pending.  This  is  an  example  of  a  successful  technology  translation  under  the  AFOSR  program. 

Even  with  improvements  in  imaging  speed,  there  is  a  drive  to  acquire  high  density 
volumetric  data  sets  in  clinical  ophthalmology  because  they  enable  the  generation  of  en  face 
OCT  images  of  the  retina.  The  acquisition  of  volumetric  data  sets  can  take  2-3  seconds.  This 
makes  volumes  susceptible  to  motion  artifacts  from  heartbeat  and  involuntary  eye  movement, 
especially  in  elderly  subjects  who  may  have  difficulties  fixating.  OCT  volumes  are  acquired  by 
raster  scanning  the  retina  where  the  OCT  beam  acquires  an  A-scan  at  one  location  on  the  retina 
at  a  time.  Transversely  scanning  the  beam  along  a  line,  called  the  fast  scan  direction,  generates  a 
B-scan  cross-sectional  image.  Acquiring  a  series  of  B-scans  at  offset  positions,  called  the  slow 
scan  direction,  generates  an  OCT  volumetric  data  set.  Motion  artifacts  occur  predominantly  in 
the  slow  scan  direction  because  the  data  acquisition  in  the  fast  scan  direction  is  rapid  enough  to 
minimize  motion  artifacts.  Working  in  collaboration  with  the  group  at  the  University  of  Erlangen, 
we  have  developed  a  software  based  registration  and  motion  correction  technique  where  two 
volumes  are  imaged  with  orthogonal  fast  scan  directions  and  then  registered  to  one  another. 
Figure  9  shows  fundus  projections  of  input  and  motion  corrected  volumes. 
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Figure  9.  Software  motion  correction.  Two  orthogonal  raster  scanned  volumes  are  acquired  and  transforms 
used  to  motion  correct  the  volumes.  The  transforms  are  iteratively  calculated  to  optimize  the  similarity 
between  the  transformed  volumes  while  penalizing  motion.  The  two  motion  corrected  volumes  are  then 
merged  to  generated  a  volume  with  improved  signal  to  noise.. 

The  motion  correction  algorithm  operated  using  two  2D  displacement  fields  for  the 
horizontal  and  vertical  raster  scanned  volumes  which  are  iteratively  optimized  to  correct  motion 
[37].  These  fields  consist  of  3D  displacement  vectors  for  every  A-scan.  A  non-linear 
optimization  algorithm  maximizes  the  similarity  between  the  displaced  volumes  while  a 
regularization  tenn  penalizes  abrupt  changes  in  the  displacement  field.  To  reduce  the  risk  of  the 
optimization  converging  at  a  local  minimum,  a  multi-resolution  approach  is  employed.  After 
registration,  the  input  volumes  are  merged  into  one  output  volume  to  improve  SNR  and  coverage 
of  the  scanned  area.  The  use  of  motion  correction  technology  to  reduce  motion  artifacts  avoids 
the  need  for  additional  hardware,  such  as  eye  tracking  devices  in  the  OCT  device.  The  motion 
correction  algorithm  can  potentially  be  applied  to  data  from  every  OCT  device  which  employs 
the  orthogonal  raster  scan  pattern. 

Recently,  an  advanced  version  of  the  motion  correction  algorithm  was  developed  [52].  The 
advanced  motion  correction  method  adds  several  elements  to  further  enhance  performance  and 
robustness  to  variations  in  input  data  quality.  Variations  in  alignment  of  the  OCT  instrument  to 
the  patient’s  eye  between  scans  can  cause  the  imaged  volumes  to  exhibit  varying  amounts  of  tilt 
along  the  transverse  directions  (Figure  10).  To  mitigate  this  effect,  the  advanced  approach 
employs  a  two-stage  optimization  (Figure  11).  The  first  stage  compensates  for  tilt  only  and  the 
adjusted  volumes  are  used  for  full  optimization  in  the  second  stage. 

Differences  in  alignment,  vignetting  or  lens  opacities  can  also  cause  differences  in  intensity 
between  subsequent  volumes  acquired  from  the  same  area.  To  correct  for  these  differences, 
illumination  correction  was  introduced  as  a  preprocessing  step.  Based  on  the  OCT  fundus 
projection  for  a  given  input  volume,  a  bias  field  is  computed  (Figure  12).  This  bias  field  is  then 
used  to  correct  the  intensities  across  the  entire  volume. 
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Figure  10.  Cross  sections  from  consecutive  orthogonal  raster  volumes  showing  different  amounts  of  tilt, 
indicated  by  the  red  lines. 


Figure  11.  Improved,  two  stage,  motion  correction  algorithm.  The  first  stage  performs  a  rough 
compensation  of  tilt  using  multi-resolution  optimization.  The  tilt-corrected  volumes  are  then  used  for  a  full 
optimization  in  stage  two. 


Input  Fundus  Bias  Field  Corrected  Fundus 


Figure  12.  An  input  OCT  fundus  image,  the  corresponding  bias  field  and  the  illumination  corrected  OCT 
fundus  image. 

Quantitative  evaluation  was  performed  using  a  large  set  of  data  acquired  from  73  healthy  and 
glaucomatous  eyes  using  SD-OCT  instruments  and  showed  that  the  advanced  motion  correction 
technology  improves  robustness  and  reproducibility  of  quantitative  measurements  between 
different  scans  when  compared  to  the  first  generation  algorithm  [52].  As  noted  previously,  IP 
was  filed  on  this  technology  and  it  was  licensed  and  transferred  to  Optovue,  Inc.  and  released  in 
a  commercial  ophthalmic  OCT  instrument  (the  Optovue,  RTVue-XR)  in  2014. 

3.  Ophthalmic  OCT 

OCT  imaging  of  the  eye  has  continued  to  be  one  of  our  program’s  major  research  areas, 
dating  back  to  our  initial  invention  and  development  of  OCT  in  1991  [53].  OCT  enables 
ophthalmologists  to  obtain  in  vivo  cross-sectional  and  volumetric  images  of  the  retina  and 
choroid  which  was  not  possible  to  obtain  with  previous  imaging  techniques.  OCT  is  now  a 
standard  imaging  modality  in  clinical  ophthalmology.  Although  there  is  considerable  commercial 
development  of  OCT  instrumentation  in  ophthalmology,  our  group  has  continued  research  in  this 
area  and  remained  at  the  forefront  of  OCT  technological  development.  In  the  previous  grant 
period,  we  developed  and  deployed  prototype,  swept  source  OCT  instruments  to  multiple  clinical 


13 


sites,  Dr.  Jay  Duker,  at  the  New  England  Eye  Center,  Dr.  Joel  Schuman  at  the  University  of 
Pittsburgh  Eye  Center,  and  Dr.  David  Huang  at  the  Oregon  Health  and  Sciences  University. 
These  instruments  used  commercially  available  swept  laser  sources  operating  at  100  kHz  A-scan 
rates  at  1050  mn  wavelength  and  were  in  use  3  years  before  the  introduction  of  commercial 
swept  source  OCT  technology. 

Working  in  collaboration  with  Praevium  Research  and  Thorlabs,  we  have  developed 
ultrahigh  speed  VCSEL  swept  light  sources  for  use  in  handheld  OCT  instrumentation,  ultrawide 
field  retinal  imaging,  and  long  range  anterior  and  full  eye  imaging.  Furthermore,  the  ultrahigh 
speed  of  the  VCSEL  light  source  is  critical  for  OCT  angiography  and  Doppler  of  the  retina. 
These  advanced  OCT  technologies  can  provide  functional  information  about  three  dimensional 
microvascular  structure  and  blood  flow  that  promise  to  elucidate  pathogenesis  as  well  as  develop 
new  diagnostic  or  progression  markers  for  ocular  diseases  which  are  leading  causes  of  blindness. 

3,1  Technology  Development  for  Clinical  Ophthalmology 

Our  group  collaborates  with  clinical  investigators  at  the  New  England  Eye  Center,  Oregon 
Health  and  Science  University  (OHSU)  and  University  of  Pittsburg  Medical  Center  (UPMC)  to 
translate  new  OCT  technology  development  in  our  laboratories  for  studies  in  the  ophthalmology 
clinic.  Technology  translation  is  a  challenging  engineering  task  because  it  requires  developing  a 
robust,  self-contained  OCT  imaging  system,  which  ergonomics  which  is  suitable  for  use  in  the 
ophthalmology  clinic  and  which  can  transported  and  operated  at  a  remote  site.  Multiple  design 
considerations  must  be  addressed  to  enable  operation  in  the  clinical  environment  which  greatly 
differs  from  the  typical  laboratory  setting.  Close  interaction  with  clinical  investigators  is 
necessary  to  understand  the  workflow  in  the  clinic,  the  clinical  applications,  and  to  design  the 
technology  so  that  it  addresses  these  needs. 

In  2011,  we  designed  and  deployed  a  1  pm  wavelength  SS-OCT  retinal  imaging  system  to 
OHSU,  shown  in  Figure  13  (A).  This  technology  employed  a  commercially  available,  short- 
cavity  swept  laser  (Axsun  Technologies,  Inc)  which  operated  at  100  kHz  A-scan  rate  [54].  The 
data  acquisition  was  optically  clocked  enabling  evenly  spaced  samples  in  wavenumber,  so  that 
the  interferometric  OCT  data  could  be  rapidly  processed  without  requiring  computationally 
costly  resampling.  The  OHSU  group,  lead  by  Dr.  David  Huang,  used  this  SS-OCT  technology 
for  their  development  of  split-spectrum  amplitude-decorrelation  angiography  (SSADA) 
algorithm  for  OCTA  [4],  This  algorithm  was  transferred  to  industry  and  has  been  integrated  into 
the  commercial  Avanti  RTVue-XR  (Optovue,  Inc.)  ophthalmic  imaging  instrument.  In  2014,  we 
designed  and  deployed  a  second  SS-OCT  instrument  to  OHSU.  This  instrument  was  especially 
designed  for  anterior  eye  imaging  and  operated  at  100  kHz  A-scan  rate  with  a  long  imaging 
range,  shown  in  Figure  13  (B).  This  system  incorporated  a  novel  optical  clock  doubling  circuit, 
similar  to  that  used  in  our  PS-OCT  studies,  to  double  the  imaging  range  of  the  system.  Doubling 
the  optical  clock  frequency  enables  both  anterior  surface  of  the  cornea  and  posterior  surface  of 
the  crystalline  lens  to  be  captured  in  a  single  image. 

In  2013,  we  designed  a  similar  100  kHz  SS-OCT  retinal  imaging  system  and  deployed  it  to 
UPMC,  Figure  13  (C).  The  UPMC  group,  led  by  Dr.  Joel  Schuman,  is  using  this  technology  for 
performing  studies  of  glaucoma.  The  combination  of  high  imaging  speed  and  the  enhanced 
penetration  of  1  pm  wavelength  enables  3D  OCT  imaging  of  the  lamina  cribrosa  structure  of  the 
optic  nerve  head.  UPMC  has  published  studies  analyzing  the  structure  of  the  lamina  cribrosa  and 
its  relation  to  glaucoma  [55,  56], 
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Our  group  has  developed  and  deployed  several  systems  to  Dr.  Duker’s  group  at  NEEC  which 
is  located  across  the  Charles  River  in  Boston.  The  system  shown  in  Figure  13  (D)  is  a  91  kHz 
ultrahigh  resolution  SD-OCT  instrument.  The  3  pm  axial  resolution  provided  by  this  prototype 
instrument  is  higher  than  the  standard  6-10  pm  in  commercial  SD-OCT  systems.  Ultrahigh 
resolution  enables  visualization  of  the  photoreceptor-RPE  complex  to  analyze  changes  in  the 
photoreceptor  outer  segments  and  RPE  which  occur  in  retinal  diseases  such  as  AMD.  We  have 
also  developed  a  100  kHz  SS-OCT  retinal  imaging  system  for  NEEC,  shown  in  Figure  13  (E), 
that  parallels  the  systems  deployed  to  the  other  collaborative  sites.  In  November  2013,  we 
developed  and  deployed  an  ultrahigh  speed  400  kHz  VCSEL  SS-OCT  retinal  imaging  system  to 
NEEC,  shown  in  Figure  13  (F).  In  addition  to  structural  imaging,  this  ultrahigh  speed 
technology  can  perform  wide  field  OCT  angiography  (OCTA)  and  acquire  multiple,  densely 
oversampled  volumes  per  cardiac  cycle  for  en  face  Doppler  imaging  of  total  retinal  blood  flow. 
These  functional  imaging  protocols  are  currently  used  to  image  patients  with  ocular  diseases 
such  as  AMD,  diabetic  retinopathy,  and  glaucoma  in  order  to  elucidate  fundamental  mechanisms 
of  disease  pathogenesis  as  well  as  to  identify  markers  for  improved  diagnosis  and  progression 
monitoring. 


Figure  13.  OCT  imaging  instruments  were  developed  and  deployed  to  the  ophthalmology  clinics  at 
multiple  collaborator  sites.  (A)  100  kHz  retinal  and  (B)  anterior  imaging  SS-OCT  systems  at  OHSU.  (C) 
100  kHz  retinal  SS-OCT  system  at  UPMC.  (D)  91  kHz  ultrahigh  resolution  retinal  SD-OCT,  (E)  100  kHz 
retinal  SS-OCT,  and  (F)  400  kHz  retinal  VCSEL  SS-OCT  at  NEEC. 
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3.2  Handheld  OCT 


Although  currently  available  commercial  OCT  systems  have  excellent  image  quality,  these 
systems  are  relatively  large  and  most  are  table  mounted,  limiting  their  use  outside  of 
ophthalmology  or  optometry  settings.  Development  of  a  handheld  OCT  technology  can  enable 
OCT  imaging  of  supine  patients  and  children,  as  well  as  facilitate  use  in  settings  such  as  primary 
care,  intraoperative,  developing  countries,  and  military  medicine.  Many  retinal  diseases  are 
asymptomatic  and  thus  go  undetected  until  irreversible  vision  loss  or  impairment  occurs  [57,  58]. 
The  rate  of  undetected  retinal  disease  is  especially  high  for  minority  communities  in  the  US  that 
do  not  have  easy  access  to  eye  care  [59,  60].  Early  detection  allows  for  early  treatment  and 
improvement  of  patient  outcome  [61,  62],  Portable  OCT  technology  would  also  be  advantageous 
for  deployment  in  remote  medical  centers,  such  as  in  field  hospitals  to  access  ocular  trauma  and 
in  third  world  countries  that  have  limited  access  to  diagnostic  technologies.  Assessment  of  OCT 
data  can  be  perfonned  remotely  by  trained  specialists  in  reading  centers. 

To  increase  the  portability  of  ophthalmic  OCT  systems,  galvanometer  scanners  used  in 
tabletop  instruments  have  been  adapted  for  handheld  instruments  in  research  [63-65]  and 
commercial  systems  such  as  the  Envisu  (Bioptigen)  and  iVue  (Optovue).  However,  OCT  beam 
scanning  required  two,  orthogonally  oriented  galvanometer  scanners  encased  in  a  bulky  metallic 
mount  for  stability  and  heat  dissipation.  To  reduce  the  size  and  weight  of  a  handheld  OCT 
instrument,  MEMs  scanning  mirrors  have  been  used  in  a  handheld  SD-OCT  system  at  70  kHz  to 
generate  2D  cross-sectional  images  for  use  in  a  primary  care  environment  [66]. 

Our  group  demonstrated  ultrahigh  speed,  swept  source  handheld  OCT  technology  using  a 
MEMS  scanning  mirror  in  combination  with  VCSEL  based  SS-OCT  [67].  The  handheld 
instrument  combined  multiple  technologies  to  overcome  many  of  the  classic  limitations  of 
handheld  imaging.  First,  we  used  a  VCSEL  light  source  to  achieve  an  imaging  speed  of  350,000 
A-scans  a  second  to  enable  volumetric  imaging  of  the  retina.  Volumetric  imaging  enables 
visualization  of  arbitrary  cross  sections  and  generation  of  en  face  images  which  can  be  used  to 
precisely  locate  retinal  features.  Second,  we  utilized  motion  correction  software  to  correct  for 
both  patient  and  operator  motion  artifacts  [37].  Third,  we  used  a  MEMS  scanning  mirror  instead 
of  galvanometers  to  scan  the  beam  and  a  plastic  3D  printed  enclosure  to  reduce  the  size  and 
weight  of  the  instrument.  Finally,  an  iris  camera,  fixation  target,  and  long  axial  imaging  range 
were  used  to  facilitate  alignment  and  reduce  imaging  session  time. 

The  system,  Figure  14,  had  an  imaging  speed  of  350,000  A-scans  per  second  and  operated  at 
1060  mn  wavelength.  The  axial  resolution  was  10  pm  in  tissue  with  a  3.1  mm  imaging  range  in 
tissue.  Two  different  handheld  enclosures  were  designed  using  3D  printing,  Figure  15.  The  first 
design  was  a  “power  grip”  style  that  has  a  handle  for  ambidextrous  operation.  The  second  design 
was  a  “camcorder”  style,  held  from  the  side  with  a  nylon  strap  supporting  the  back  of  the  hand. 
A  rotational  adjustment  was  used  to  translate  the  ocular  lens  to  compensate  refractive  errors  of  - 
12  to  +6  Diopters.  Both  enclosures  had  identical  optical  components,  including  an  iris  camera 
and  manually  adjustable  fixation  target.  The  iris  camera  provided  a  video  of  the  pupil  and  iris  to 
aid  in  aligning  the  OCT  beam  through  the  pupil.  The  manual  fixation  target  could  move 
horizontally  for  OCT  imaging  of  the  macula  or  optic  nerve  head.  The  video  from  the  iris  camera 
and  preview  OCT  scans  were  displayed  on  a  3.5  inch  LCD  screen  on  the  device. 
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Figure  14.  Swept  Source  OCT  imaging  engine  using  a  high  speed  VCSEL  light  source  at  1060  nm.  A 
calibration  MZI  is  used  to  generate  an  optical  clock  signal.  MEMS  scanning  mirror  (SM).  Dichroic  mirror 
(DM).  Polarization  controller  (PC).  Dispersion  compensation  glass  (DC). 


Figure  15.  Photographs  of  the  (A,  B)  power  grip  style  and  (C,  D)  camcorder  style  designs. 
The  optical  components  inside  both  enclosures  are  identical. 
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Figure  16.  Internal  optical  layout  of  handheld  OCT  and  unfolded  optical  design  showing  the  (A)  OCT 
1060  nm  optical  path,  (B)  iris  camera  optical  path,  and  (C)  fixation  target  optical  path. 
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Figure  16  shows  the  OCT  (A),  iris  camera  (B),  and  fixation  target  (C)  optical  pathways.  A 
2.4  mm  diameter  2D  gimbal-less  tip-tilt  MEMS  mirror  (Mirrorcle  Technologies)  was  used  to 
scan  the  OCT  beam.  Due  to  the  low  resonance  frequency  of  the  MEMS  mirror,  volumetric 
images  were  acquired  by  scanning  the  mirror  sinusoidally  in  one  direction  and  linearly  in  the 
other  direction.  Two  orthogonally  scanned  volumes  were  acquired,  motion  corrected  and  merged 
using  our  software  motion  correction.  Figure  17  shows  a  10  mm  x  10  mm  volumetric  image 
taken  with  the  handheld  system.  Two  1400  x  350  A-scan  sinusoidal  raster  volumes  were 
acquired  in  1.4  seconds  per  volume.  The  two  volumes  were  converted  to  350  x  350  A-scan  linear 
volumes  and  motion  corrected  to  generate  a  merged  350  x  350  A-scan  volume.  This  volume  was 
centered  between  the  fovea  and  optic  nerve  head  by  adjusting  the  fixation  target.  Since  this  wide 
field  volumetric  scan  captures  the  macular  and  optic  nerve  head  area,  the  operator  does  not  need 
to  image  both  areas  separately. 


Figure  17.  Hand  held  OCT  imaging.  Motion-corrected,  wide  field  10  mm  x  10  mm,  350  x  350  A-scan 
volume  generated  from  two  volumes  acquired  in  1.4  seconds  each.  (A)  En  face  OCT  fundus  image.  (B) 
Original  vertical  cross-sectional  image  from  the  horizontal  raster  scanned  volume.  (C)  Vertical  cross- 
sectional  image  from  the  motion-corrected  merged  volume.  (C,  D)  Colors  indicate  cross  section  locations 
on  the  fundus.  (E)  Interpolated  3.4  mm  diameter  circumpapillary  scan  extracted  from  motion-corrected 
volumetric  data.  Scale  bars:  1  mm. 


3.3  Wide  Field  Retinal  Imaging 

For  volumetric  OCT  imaging,  the  A-scan  sampling  interval  should  be  smaller  than  the 
illumination  spot  diameter  in  order  to  ensure  optimal  resolution  in  order  to  visualize  focal 
pathology.  Therefore,  faster  speeds  are  required  to  image  wider  areas,  since  the  data  acquisition 
times  for  retinal  imaging  are  limited.  The  400  kHz  ultrahigh  speed  VCSEL  swept  source  OCT 
technology  enables  much  wider  retinal  imaging  fields,  while  still  maintaining  comprehensive 
structural  information  [4].  Figure  18  shows  a  14  mm  x  14  mm,  800  x  800  A-scan  motion- 
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corrected  volume  of  the  retina  generated  from  two  orthogonally  scanned  volumes  each  acquired 
in  ~3  seconds.  Wide  field  imaging  enables  coverage  of  both  the  macula  and  the  optic  nerve  head 
in  a  single  data  acquisition,  simplifying  workflow  and  reducing  the  time  required  to  image 
patients. 


Figure  18.  Wide  field  volumetric  OCT  imaging  of  the  human  eye.  (A)  OCT  fundus  image  from  a  large  14 
mm  x  14  mm,  800^800  A-scan  motion-corrected  volumetric  data  set.  (B-D)  Arbitrary  cross-sectional 
images  can  be  extracted  from  volumetric  data.  The  data  contains  comprehensive  structural  information  on 
both  retina  and  choroid.  Scale  bar:  1  mm. 


3.4  Long  Range  Imaging  of  the  Anterior  Eye  and  Full  Eye  Length  Imaging 

Utilizing  the  long  imaging  range  and  adjustability  of  the  MEMS-tunable  VCSEL,  we 
developed  prototype  multi-modal  SS-OCT  technology  that  can  be  adapted  for  visualizing  the 
human  retina,  the  anterior  segment  of  the  eye,  and  the  full  eye  length,  with  minimal  hardware 
modification  in  the  patient  interface  [4],  For  retinal  imaging,  the  A-scan  rate  was  set  to  580  kHz 
with  1.9  mm  imaging  range  in  tissue.  For  imaging  the  anterior  segment,  the  A-scan  rate  was  set 
to  100  kHz  with  13.6  mm  imaging  range.  In  order  to  achieve  38  mm  imaging  range  to  cover  the 
full  eye  length,  the  A-scan  rate  was  further  reduced  to  50  kHz.  Figure  19  shows  a  3D  rendering 
and  representative  cross  sectional  image  from  an  OCT  volume  of  the  entire  anterior  segment 
from  the  anterior  cornea  to  the  posterior  surface  of  the  ocular  lens.  Previously,  Visualizing  the 
entire  anterior  segment  with  OCT  has  been  challenging  because  other  SS-OCT  light  sources 
have  limited  coherence  lengths  and  limited  imaging  range. 

The  measurement  of  the  optical  surfaces  of  the  anterior  eye  is  essential  for  monitoring 
refractive  or  cataract  surgery  procedures  such  as  laser-assisted  in  situ  keratomileusis  (LASIK), 
photorefractive  keratocomy  (PRK),  phototheraputic  keratocomy  (PTK),  intracorneal  rings, 
astigmatic  keratocomy,  lamellar  keratoplasty  and  intraocular  lens  implants  [68-71].  The  ultrahigh 
imaging  speed  and  long  imaging  range  of  VCSEL  SS-OCT  in  combination  with  software  motion 
correction  enable  in  vivo  anterior  segment  biometry.  Figure  20  shows  examples  of  biometric 
measurements  of  the  corneal  and  the  crystalline  lens  surfaces.  Image  distortion  caused  by  light 
refraction  was  corrected  using  a  3D  ray-tracing  algorithm.  Advanced  motion  correction  methods 
are  crucial  for  performing  keratometry  with  OCT  because  mapping  optical  power  requires 
accurate  measurement  of  the  second  derivative  of  the  surface  contour. 
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The  unprecedented  imaging  range  of  the  VCSEL  also  enables  full-eye  OCT,  as  shown  in 
Figure  21.  Full-eye  OCT  spans  the  full  eye  length  from  the  cornea  to  the  retina,  enabling  precise 
measurement  of  intraocular  distances  [72].  We  have  investigated  the  measurement 
reproducibility  of  the  full-eye  OCT  ocular  biometry,  demonstrating  good  reproducibility  and 
correlation  with  commercially  available  optical  and  ultrasonic  biometers  [51]. 


Figure  19.  Anterior  segment  imaging  with  high  speed,  long  range  SS-OCT.  (A)  Motion-corrected  volume 
comprising  400x400  A-scans  over  a  15  mm  x  15  mm  area.  (B)  Representative  central  cross-sectional 
image  from  the  volumetric  data  set. 
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Figure  20.  Examples  of  anterior  eye  biometry.  Top  row  shows  the  topographic  map  of  the  anterior 
chamber  surfaces  and  the  bottom  row  shows  the  residual  elevation  from  spherical  fit. 


Figure  21.  Full-eye  OCT  can  visualize  the  anterior  and  posterior  eye  in  a  single  volume.  (A)  Full-eye  OCT 
volume  with  500  x  500  A-scans,  15  mm  x  15  mm  area  and  42  mm  imaging  range.  (B)  Central  cross  section 
and  averaged  intensity  profile  of  extracted  central  A-scans  shows  the  axial  distances  of  ocular  surfaces. 
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3.5  OCT  Angiography  of  the  Retina  and  Choroid 

Current  Techniques  for  Visualizing  Retinal  and  Choriocapillaris  Microvasculature 

The  retina  and  choroid,  which  are  the  two  blood  supplies  to  the  eye,  are  thought  to  play  an 
important  role  in  a  variety  of  ocular  diseases,  including  age-related  macular  degeneration  (AMD) 
and  diabetic  retinopathy.  Traditionally,  ocular  vasculature  has  been  visualized,  both  in  the 
research  setting  and  clinically,  using  fluorescein  angiography  (FA)  and  indocyanine  green 
angiography  (ICGA)  which  require  intravenous  dye  injection.  While  FA  can  visualize  retinal 
vasculature,  visualization  of  the  choriocapillaris/choroid  using  FA  is  challenging  because  the 
blue-green  excitation  wavelength  of  fluorescein  is  partially  absorbed  by  the  macular  xanthophyll 
and  retinal  pigment  epithelium  (RPE),  and  -20%  of  the  injected  fluorescein  does  not  bind  to 
albumin  and  leak  outs  of  the  choriocapillaris  fenestrations,  creating  diffuse  hyperfluorescence 
which  obscures  vasculature  [73].  In  contrast,  indocyanine  green  angiography  (ICGA)  enables 
visualization  of  choroidal  circulation  because  its  near  infrared  excitation  wavelengths  are  not  as 
readily  absorbed  by  the  macular  xanthophyll  and  RPE,  and  -98%  of  the  injected  indocyanine 
green  is  bound  to  plasma  protein,  preventing  it  from  leaking  out  of  the  choriocapillaris 
fenestrations  [73].  ICGA  has  also  been  shown  to  be  capable  of  visualizing  the  choriocapillaris 
circulation  [74].  However,  since  ICGA  is  not  depth  resolved,  the  task  of  separating 
choriocapillaris  blood  flow  from  that  of  deeper  choroidal  vasculature  is  complicated  and  may 
only  be  possible  if  it  is  assumed  that  the  blood  flow  velocity  of  the  larger  choroidal  arteries  is 
small  compared  to  the  choriocapillaris  [74,  75].  In  addition  to  the  limited  ability  of  FA  and  ICGA 
to  image  the  choriocapillaris,  both  require  the  injection  of  exogenous  dyes,  making  the  imaging 
procedure  invasive  and  time  consuming. 

OCT  angiography  (OCTA)  can  generate  three-dimensional  microvascular  angiograms  in 
vivo,  rapidly  and  without  requiring  injection  of  exogenous  dyes  [7-14,  76].  OCTA  is  based  on  the 
fact  that  stationary  tissue  generates  a  time-independent  B-scan  image,  while  flowing  blood 
motion  produces  time-dependence  in  the  B-scans.  If  repeated  B-scans  of  stationary  tissue  are 
acquired  at  the  same  location  over  time,  then  each  B-scan  will  be  identical.  However,  if  there  is 
motion,  for  example  flowing  erythrocytes,  then  repeated  B-scans  will  differ.  These  differences 
can  be  quantitatively  represented  by  a  decorrelation  signal,  where  large  differences  in  image 
features  (corresponding  to  fast  flow)  generate  high  decorrelation  and  small  differences 
(corresponding  to  slow  flow)  generate  low  decorrelation.  A  number  of  different  techniques  for 
quantitating  the  decorrelation  signal  have  been  demonstrated  [5,  10,  77].  Since  OCT  angiograms 
are  obtained  from  repeated  intensity  B-scans,  they  are  intrinsically  co-registered  with  structural 
data.  OCTA  therefore  enables  simultaneous  assessment  of  structure  and  micro  vasculature. 

Compared  with  conventional  dye-based  angiography  techniques,  OCTA  offers  several 
important  advantages.  First,  OCTA  does  not  require  the  injection  of  exogenous  dyes  and  is 
noninvasive.  Second,  unlike  dye-based  angiography,  OCTA  is  fast  (typical  acquisition  time  <  4 
seconds),  and  can  be  performed  repeatedly,  potentially  during  every  patient  visit.  Third,  OCTA 
enables  depth  resolved  imaging  of  both  the  retinal  and  choroidal  vasculature.  Taken  together, 
these  advantages  make  OCTA  a  powerful  technique  for  studying  retinal  and  choroidal 
microvasculature,  particularly  in  longitudinal  clinical  studies. 

Development  of  SS-OCT  Angiography  using  the  VCSEL  Swept  Light  Source 
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Our  group  has  demonstrated  ultrahigh  speed  swept  source  OCT  (SS-OCT)  for  ophthalmic 
applications  using  VCSEL  technology  and  developed  a  phase-stabilized  prototype  system  with  a 
400  kHz  A-scan  rate,  approximately  5-10  times  faster  than  standard  commercial  systems  [5],  The 
instrument  used  a  VCSEL  centered  at  -1050  nm  wavelength  with  a  -80  nm  tuning  range.  OCT 
signals  are  sampled  with  a  high-speed  digitizer  externally  clocked  at  up  to  — 1.1  GHz  using  a 
Mach-Zehnder  interferometer  calibration  signal  yielding  an  imaging  range  of  -2.1  mm  in  tissue. 
Imaging  at  -1  pm  wavelengths  has  the  advantage  over  840  nm  wavelength  light  sources  in 
current  commercial  SD-OCT  systems  in  that  it  is  less  sensitive  to  attenuation  by  ocular  opacities 
and  allows  deeper  light  penetration  into  tissue,  such  as  the  RPE  and  choroid.  The  ultrahigh  speed 
is  especially  advantageous  for  imaging  the  choriocapillaris,  which  has  very  features  requiring  a 
high  A-scan  sampling  density. 

Demonstration  of  SS-OCT  Angiography  Imaging  in  Normal  Subjects 

Our  group  demonstrated  in  vivo  OCTA  imaging  of  the  choriocapillaris  and  choroidal 
microvasculature  in  a  cohort  of  normal  human  subjects  using  OCTA  [78].  The  vascular  patterns 
of  the  OCTA  were  shown  to  be  consistent  with  those  seen  in  previous  histological  and  electron 
micrograph  studies  (Figure  22)  [17,  79]. 
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Figure  22.  1.5  mm  x  1.5  mm  en  face  OCT  angiograms  at  three  different  fundus  locations.  (A,  D,  G)  OCT 
angiogram  of  the  choriocapillaris,  OCT  angiogram  of  the  Sattler’s  layer,  and  en  face  OCT  intensity  image 
(slice)  corresponding  to  the  choriocapillaris  layer,  respectively,  at  the  central  fovea.  (B,  E,  H)  OCT 
angiogram  of  the  choriocapillaris,  OCT  angiogram  of  the  Sattler’s  layer,  and  en  face  OCT  intensity  image 
(slice)  corresponding  to  the  choriocapillaris  layer,  respectively,  at  6  mm  temporal  to  the  central  fovea.  The 
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arrowhead  indicates  venules  and  arrow  arteriole  in  (E).  (C,  F,  I)  OCT  angiogram  of  the  choriocapillaris, 

OCT  angiogram  of  the  Sattler’s  layer,  and  en  face  OCT  intensity  image  (slice)  corresponding  to  the 
choriocapillaris  layer,  respectively,  at  12  mm  nasal  to  the  central  fovea.  (J,  K)  Electron  micrographs  of 
corrosion  casts  reproduced  from  Olver  et  al.  with  permission.  Scale  bars:  250  mm. 

Exudative  AMD 

Exudative  age-related  macular  degeneration  (AMD),  a  pathology  characterized  by  choroidal 
neovascularization  (CNV),  is  a  leading  cause  of  vision  loss  and  impairment  in  developed 
countries.  Working  in  collaboration  with  Dr.  Jay  Duker  at  the  NEEC,  we  performed  a  clinical 
study  of  OCTA  imaging  and  analysis  of  17  eyes  with  exudative  AMD  with  active  CNV  (15 
patients;  79.7  ±  8.3  years  old).  Enrollment  is  ongoing.  Choroidal  neovascularization  was  clearly 
visualized  in  16  of  the  19  eyes  with  exudative  AMD.  In  all  of  the  16  eyes  in  which  CNV  was 
visualized,  the  CNV  was  located  above  regions  of  severe  choriocapillaris  alteration.  In  14  of 
these  eyes,  the  CNV  lesions  were  surrounded  by  regions  of  choriocapillaris  alteration.  Figure  23 
shows  OCT  and  OCTA  images  from  an  87  year  old  treatment  naive  exudative  AMD  patient  with 
a  minimally  classic  CNV.  All  images  were  obtained  from  a  single  volumetric  scan.  The  results  of 
this  study  suggest  mechanisms  for  CNV  development  and  its  link  to  choriocapillaris  alterations. 
The  study  also  suggests  a  role  for  OCTA  in  the  identification  of  CNV  for  determining  treatment 
of  exudative  AMD.  The  results  of  our  exudative  AMD  studies  were  recently  featured  on  the 
cover  of  Ophthalmic  Surgery,  Lasers  and  Imaging  Retina,  Nov/Dec  2014  [1], 

Dry  AMD 

While  the  most  severe  vision  loss  has  historically  been  associated  with  the  exudative  or  wet 
form  of  AMD  which  is  characterized  by  the  formation  of  choroidal  neovascularization  (CNV), 
the  atrophic  or  dry  form  of  AMD  which  is  characterized  by  geographic  atrophy  (GA),  is  likely  to 
become  the  most  common  cause  of  severe  loss  in  the  future.  This  shift  is  largely  due  to  the 
success  of  vascular  endothelial  growth  factor  inhibitors.  We  have  performed  a  clinical  study  of 
OCTA  imaging  and  analysis  of  12  eyes  with  non-exudative  AMD  with  GA  (7  patients;  75.9  ± 
6.1  years  old).  Enrollment  is  ongoing.  In  12/12  eyes  OCTA  revealed  evidence  of  CC  atrophy 
within  the  regions  of  GA.  In  10/12  eyes  OCTA  showed  CC  flow  impairment  extending  beyond 
the  margin  of  GA.  In  14/16  eyes  choriocapillaris  alteration  was  observed  in  regions  extending 
beyond  the  CNV  lesion.  Figure  24  illustrates  the  ability  of  OCTA  to  examine  the 
choriocapillaris  in  a  75  year  old  atrophic  age-related  macular  degeneration  (AMD)  patient.  These 
studies  also  investigated  the  role  of  OCTA  imaging  parameters  in  detennining  the  sensitivity  of 
OCTA  to  low  flows  and  provided  a  method  to  differentiate  choriocapillaris  flow  impairment 
from  atrophy.  More  extensive  cross-sectional  studies  are  ongoing  and  are  needed  to  investigate 
the  potential  of  OCTA  for  monitoring  disease  progression  and  response  to  therapy.  Recruitment 
of  age  matched  normal  subjects  will  be  especially  important  in  order  to  separate  effects  of 
nonnal  aging  from  alterations  due  to  disease. 
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Figure  23.  An  87  year  old  patient  with  a  CNV  lesion  exhibiting  both  occult  and  classic  components.  (A) 
OCT  intensity  projection.  (B)  En  face  OCTA  projection  showing  retinal  vasculature.  (C)  En  face  slice  of 
the  OCT  intensity  volume  at  the  depth  of  the  choriocapillaris.  The  white  dashed  contour  shows  an  area  of 
lower  intensity.  (D)  En  face  slice  of  the  OCTA  volume  at  the  depth  of  the  choriocapillaris.  (E)  En  face 
projection  of  the  OCTA  volume  through  the  depths  spanned  by  the  CNV.  The  red  shading  corresponds  to 
the  occult  components  of  the  lesion  while  the  yellow  shading  corresponds  to  the  classic  components.  The 
arrows  point  to  motion  artifacts  that  appear  as  horizontal  lines  through  the  image.  (F)  Early  phase  FA. 
Arrows  point  to  the  classic  component.  (G)  Early  phase  ICGA.  Arrows  point  to  the  occult  component.  (H) 
Late  phase  FA.  Arrows  point  to  area  of  stippled  hyperfluorescence.  (I)  OCT  retinal  thickness  map.  (J,  K) 
OCT  intensity  and  OCTA  cross-sections,  respectively,  corresponding  to  the  white  dashed  horizontal  lines 
in  (E).  The  bracket  in  (J)  roughly  spans  the  depths  in  which  the  OCTA  B-scan  is  uninterpretable  due  to 
shadowing  from  the  choriocapillaris/choroidal  vasculature.  The  yellow  arrows  in  (K)  indicate  the  choroidal 
thickness  away  from  the  lesion;  the  white  arrows  show  an  area  of  shadowing;  and  the  blue  arrows 
correspond  to  a  region  of  choroidal  atrophy.  All  OCT  images  are  from  a  6  mm  x  6  mm  area  and  all  scale 
bars  are  1  mm.  All  OCT  and  OCTA  images  were  generated  from  a  single  volumetric  scan. 


24 


Figure  24.  Fundus  autofluorescence  (FAF)  image,  OCT  and  OCTA  in  a  75  year  old  atrophic  age-related 
macular  degeneration  (AMD)  patient  with  a  visual  acuity  of  20/20.  The  FAF  (A)  and  the  en  face  OCT  (B) 
clearly  show  the  region  of  GA,  outlined  by  the  yellow  dashed  line  in  the  en  face  OCT.  The  GA  region 
appears  lighter  due  to  increased  light  penetration  into  the  choroid  caused  by  RPE  atrophy.  The  white  arrow 
shows  the  region  of  foveal  sparing.  (C)  An  OCTA  of  the  retinal  microvasculature,  generated  by  mean 
projection  of  the  OCTA  volume  between  the  inner  limiting  membrane  (ILM)  and  retinal  pigment 
epithelium  (RPE).  The  retinal  vasculature  appears  normal.  (D)  An  en  face  OCTA  of  a  4.4  pm  thick  slab  at 
the  level  of  the  choriocapillaris  (CC)  using  a  ~7.5  ms  interescan  time.  The  yellow  dashed  contour  from  (B) 
is  superposed,  and  severe  CC  alteration  appears  within  this  contour.  Flow  in  the  area  of  foveal  sparing, 
indicated  by  the  yellow  arrow,  is  also  visible.  Severe  CC  alteration  is  evident  outside  the  GA  margin.  (E) 
An  en  face  OCT  of  the  same  4.4  pm  thick  slab  at  the  level  of  the  CC  as  in  (D),  using  a  ~15  ms  interscan 
time  which  has  increased  sensitivity.  Magnified  views  of  the  solid  orange  and  green  boxes  of  (D)  and  (E) 
are  shown  in  (F)  and  (G),  respectively.  Note  that  some  choroidal  vessels  that  are  not  visible  in  (F)  become 
visible  in  the  high  sensivity  image  (G).  Magnified  views  of  the  dashed  orange  and  green  boxes  of  (D)  and 
(E)  are  shown  in  (FI)  and  (I).  OCT  intensity  (top)  and  OCTA  (bottom)  B-scans  through  the  red,  blue  and 
purple  horizontal  dashed  lines  in  (D)  are  shown  in  (J),  (K)  and  (L),  respectively.  All  scale  bars  are  1  mm. 


Diabetic  Retinopathy 

Both  retinal  and  choroidal  blood  circulation  are  believed  to  be  altered  in  diabetic  retinopathy 
(DR),  a  leading  cause  of  vision  loss  or  impairment  in  developed  countries.  There  are  multiple 
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vascular  abnormalities  that  are  potential  early  markers  for  DR.  Early  changes  in  non-proliferative 
diabetic  retinopathy  (NPDR)  include  microaneurysms,  foveal  avascular  zone  (FAZ) 
enlargement,  and  capillary  dropout.  We  have  performed  a  clinical  study  of  OCTA  imaging  and 
analysis  of  9  eyes  with  proliferative  diabetic  retinopathy  (PDR),  29  eyes  with  non-proliferative 
diabetic  retinopathy  (NPDR),  and  53  diabetic  eyes  without  retinopathy.  Enrollment  is  ongoing. 
Figure  25  illustrates  the  ability  of  OCTA  to  visualize  microvascular  changes  of  a  55  year  old 
diabetic  patient  with  PDR.  Preliminary  results  suggest  that  microvascular  changes  appear  in 
patients  without  overt  retinopathy. 


Figure  25.  (A)  FA,  (B)  fundus  photograph,  (C)  enlarged  view  of  the  FA,  and  (D)  OCT  angiogram  from  a 
55  year  old  diabetic  patient  with  PDR.  The  enlarged  area  indicated  by  the  orange  square  in  (A)  is  shown  in 
(C).  The  OCT  angiogram  in  (D)  show  high  correspondence  to  the  FA  in  (C).  Many  of  the  microaneurysms 
indicated  by  the  yellow  circles  appear  as  hyperfluorescent  spots  in  the  FA  while  they  appear  as  abnormal 
capillary  loops  in  the  OCT  angiogram.  OCT  angiogram  acquired  from  a  3  mm  x  3  mm  area. 

3.6  En  Face  Doppler  OCT  Measurement  of  Total  Retinal  Blood  Flow 

Measurement  of  Total  Retinal  Blood  Flow 

Measurement  of  retinal  blood  flow  in  vivo  is  an  important  research  area  in  ophthalmic 
imaging  since  alterations  in  retinal  perfusion  are  suggested  to  play  an  important  role  in 
pathogenesis  of  ocular  diseases  such  as  glaucoma  [80-83]  and  diabetic  retinopathy  [84-87]. 
Quantitative  assessment  of  total  retinal  blood  flow  (TRBF)  enables  objective  comparison  of 
retinal  circulation  among  different  subjects  and  simplifies  longitudinal  tracking  of  changes 
associated  with  disease  progression.  TRBF  measurement  can  provide  useful  information  for 
investigating  ocular  pathophysiology  as  well  as  for  identifying  potential  early  diagnostic 
markers. 

Conventional  Techniques  of  TRBF  Measurement 

Multiple  imaging  modalities  such  as  video  fluorescein  angiography,  laser  speckle  flowgraphy, 
laser  Doppler  velocimetry,  scanning  laser  Doppler  flowmetry,  and  ultrasound  color  Doppler 
imaging  can  be  used  to  investigate  quantitative  retinal  blood  flow  [88-94].  However,  objective 
measurement  of  TRBF  remains  challenging.  Video  fluorescein  angiography  and  laser  speckle 
flowgraphy  cannot  directly  measure  blood  flow  velocity,  but  yield  indirect  metrics  such  as 
arteriovenous  passage  time  and  mean  blur  rate.  Scanning  laser  Doppler  flowmetry  (e.g. 
Heidelberg  Retina  Flowmeter,  Heidelberg  Engineering,  Heidelberg,  Germany)  can  map  the  flow 
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velocity  component  along  the  probe  beam  direction  over  a  two-dimensional  scan  area;  however, 
the  absolute  flow  velocity  is  required  for  quantitative  flow  rate  measurement.  Ultrasound  color 
Doppler  imaging  can  measure  flow  velocity  at  different  tissue  depths  but  has  a  ~200  um  spatial 
resolution  which  is  insufficient  to  visualize  detailed  structure  of  the  central  retinal  vasculature, 
complicating  accurate  blood  flow  calculation. 

Among  these  techniques,  bidirectional  laser  Doppler  velocimetry  (BLDV)  combined  with 
fundus  photography  is  one  of  the  most  well  suited  methods  for  quantitative  TRBF  measurement. 
This  technique  measures  maximum  absolute  flow  velocity  in  a  vessel  by  detecting  the  Doppler 
frequency  shifts  of  the  light  backscattered  in  two  different  directions.  Blood  flow  in  the  vessel  is 
calculated  as  the  product  of  the  maximum  absolute  flow  velocity  and  the  effective  vessel  cross- 
sectional  area,  estimated  from  the  vessel  diameter  measured  in  a  fundus  photograph.  TRBF  is 
then  calculated  as  the  sum  of  the  blood  flow  in  major  retinal  vessels  near  the  optic  disc.  BLDV 
measurement  of  TRBF  was  validated  in  healthy  subjects  [95-98]  and  used  for  investigating 
TRBF  in  diabetic  patients  [99-102].  However,  the  calculation  of  effective  vessel  cross-section 
area  relies  on  the  assumption  that  vessel  shapes  are  round  and  blood  flow  velocity  follows  a 
parabolic  profile.  This  assumption  can  be  inaccurate  in  some  subjects  and  may  result  in  flow 
measurement  errors.  Moreover,  measurement  times  can  be  preclusive  because  individual  major 
retinal  vessels  must  be  measured  separately,  requiring  20  to  55  minutes  per  eye  [103].  Due  to  this 
limitation,  the  number  of  patients  in  most  TRBF  studies  using  BLDV  has  been  limited. 

En  face  Doppler  OCT 

Fourier-domain  OCT  has  enabled  phase-sensitive  imaging  methods  such  as  Doppler  OCT, 
which  promises  simpler  and  more  robust  measurements  of  TRBF  using  depth-resolved  vessel 
structure  and  Doppler  velocity  [104].  Conventional  Doppler  OCT  techniques  require  measuring 
the  angle  of  blood  vessel,  because  Doppler  OCT  detects  only  the  axial  velocity  component. 
Previous  Doppler  OCT  methods  have  used  double  circumpapillary  scans  that  intercept  central 
retinal  vasculature  to  detect  the  axial  velocity  and  measure  the  Doppler  angle,  the  angle  between 
the  OCT  beam  and  blood  flow  direction  for  each  major  retinal  vessel  [105,  106].  However,  this 
technique  was  sensitive  to  small  errors  in  Doppler  angles,  making  automatic  measurements 
challenging.  Our  collaborators  from  Dr.  David  Boas’  group  at  MGH  developed  a  new  en  face 
Doppler  OCT  technique  for  measuring  blood  flow  in  the  small  animal  brain,  which  does  not 
require  angle  measurements  [18].  En  face  Doppler  measures  blood  flow  by  volumetric  raster 
scanning  the  vessels  and  integrating  the  axial  velocity  components  in  an  en  face  plane.  The  en 
face  velocity  can  be  integrated  to  obtain  the  total  flow  without  requiring  the  Doppler  angle.  Our 
group  used  a  200  kHz  SS-OCT  instrument  to  demonstrate  that  human  total  retinal  blood  flow  can 
be  measured  by  volumetrically  scanning  a  small  area  at  the  optic  disc  [107].  This  is  an  example 
of  translating  imaging  techniques  developed  in  the  brain  to  the  retina  and  shows  the  synergy 
between  methods  in  brain  science  and  ophthalmology. 

En  face  Doppler  OCT  with  High-Speed  OCT  Technology 

Our  group  continued  development  en  face  Doppler  OCT  measurement  techniques  using 
ultrahigh  speed  SD-OCT  and  SS-OCT  to  resolve  pulsatile  retinal  blood  flow  and  improve 
measurement  reproducibility  in  small  animals  and  human  subjects  [5,  108].  Small  animal  studies 
are  particularly  challenging  because  heart  rates  in  rats  can  be  300  beats  per  minute.  Ultrahigh 
speed  en  face  Doppler  OCT  enabled  volume  acquisition  rates  of  55  Hz  in  small  animals  using 
SD-OCT  and  7.6  Hz  in  humans  using  SS-OCT.  These  volume  acquisition  rates  are  several  times 
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faster  than  the  normal  heart  rates  in  small  animals  and  humans,  enabling  the  acquisition  of 
multiple  flow  measurements  per  cardiac  cycle.  Figures  26  (B-F)  show  an  example  of  pulsatile 
retinal  blood  flow  measurement  in  a  normal  human  subject  performed  with  an  ultrahigh  speed 
swept  source  OCT  prototype  using  a  1  pm  wavelength  VCSEL  light  source  at  a  400  kHz  A-scan 
rate.  Each  flow  value  measurement  in  Figure  26  (B)  requires  a  3D  OCT  volume. 


Figure  26.  En  face  Doppler  OCT  provides  a  robust  and  fully  automatable  measurement  of  total  retinal 
blood  flow.  (A)  Flow  is  measured  in  en  face  plane  avoiding  the  need  to  measure  the  Doppler  angle  a.  (B) 
Ultrahigh  speed  OCT  at  1050  nm  enables  repeated  volumetric  imaging  (750  x  60  A-scan  volumes  at  7.6  Hz 
rate)  to  measure  pulsatile  flow.  (C-F)  Pulsatile  flow  in  en  face  plane  of  central  retinal  vasculature.  Scale  bar: 

500  pm. 

Synchronized  En  face  Doppler  OCT  with  Pulse  Oximetry  Cardiac  Gating 

In  order  to  enable  larger  scale  clinical  studies,  we  developed  a  cardiac  gating  method  for 
enhanced  sampling  of  blood  flow  pulsatility  in  en  face  Doppler  with  commercially  available  100 
kHz  A-scan  rate  SS-OCT  technology.  A  pulse  oximetry  plethysmogram  was  obtained  from 
subject’s  earlobes  and  was  connected  to  the  OCT  system  to  synchronize  the  start  of  a  multiple 
volume  OCT  acquisition  sequence  with  the  subject’s  cardiac  cycle.  Total  flow  was  measured 
using  three  synchronized,  multiple  volume  acquisitions.  Highly  repeatable  mean  total  retinal 
blood  flow  measurement  with  fully  automatic  processing  was  made  possible  using  specialized 
segmentation  algorithm  and  retrospective  correction  of  cardiac  gating.  In  a  repeatability  study  of 
mean  total  retinal  blood  flow  measurement  on  10  healthy  young  subjects,  we  observed  an 
improvement  from  8.0%  median  coefficient  of  variance  (COV)  to  4.9%  COV  using  cardiac 
gating.  In  the  same  study,  ultrahigh  speed  asynchronous  acquisition  at  400  kHz  A-scan  rate  and 
7.6  Hz  volume  rate  yielded  median  COV  of  6.1%.  These  results  demonstrate  that  the  simple 
addition  of  pulse  oximeter  can  enable  repeatable  measurement  of  total  retinal  blood  flow  in  the 
clinic  using  commercially  available  OCT  imaging  speeds. 

Clinical  Studies  of  Total  Retinal  Blood  Flow 
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In  an  ongoing  clinical  study  at  NEEC  using  the  400  kHz  prototype  SS-OCT  instrument,  we 
measured  total  retinal  blood  flow  in  18  eyes  with  non-proliferative  diabetic  retinopathy,  4  eyes 
with  proliferative  diabetic  retinopathy,  16  diabetics  without  clinically  evident  retinopathy,  21 
eyes  with  primary  open-angle  glaucoma,  and  3  eyes  with  normal  tension  glaucoma.  The  study  is 
ongoing,  however  preliminary  results  suggest  that  patients  with  primary  open-angle  glaucoma 
exhibit  lower  mean  retinal  blood  flow  than  normal  subjects.  The  mean  flow  difference  between 
nonnal  eyes  and  eyes  with  primary  open-angle  glaucoma  was  not  large  compared  to  the  variation 
inside  each  group,  suggesting  that  the  statistical  significance  of  the  mean  flow  difference  has  to 
be  further  tested.  In  eyes  with  diabetic  retinopathy,  blood  flow  had  a  wide  range  of  variation, 
occasionally  exhibiting  both  lower  [92-94]  and  higher  [87,  99,  100]  flows  than  the  normal  range. 
Previous  studies  in  diabetic  retinopathy  have  also  reported  both  an  increase  and  decrease  in 
blood  flow,  depending  on  the  severity  of  the  disease.  Therefore,  further  investigations  involving 
detailed  classification  of  retinopathy  severity  and  its  correlation  with  blood  flow  alteration  are 
important. 

4.  Brain  Imaging 

4,1  Introduction  to  OCT  for  Neuroscience  Applications 

In  collaboration  with  Dr.  David  Boas  at  Martinos  Center  of  Biomedical  Imaging, 
Massachusetts  General  Hospital,  Harvard  Medical  School,  we  have  developed  advanced  OCT 
technologies  for  neuroscience  applications.  Many  of  these  methods  can  be  translated  to 
fundamental  and  clinical  research  in  ophthalmology.  These  studies  are  also  relevant  for 
investigating  the  pathogenesis  of  traumatic  brain  injury  (TBI),  a  problem  of  direct  military 
relevance.  Dr  Boas’s  group  has  extensive  infrastructure  for  Magnetic  Resonance  Imaging  as  well 
as  OCT  and  a  track  record  of  collaborations  with  neurologists  specializing  in  brain  science, 
neurovascular  coupling,  brain  trauma,  stroke  and  migraine. 

Optical  imaging  methods  have  had  a  significant  impact  on  the  field  of  neuroimaging  over  the 
past  20  years.  Optical  imaging  is  now  widely  used  in  studies  of  cerebrovascular  and  cellular 
physiology  and  pathology.  While  macroscopic  methods  using  diffuse  light  (optical  intrinsic 
signal  imaging  [109],  laser  Doppler  imaging  [110],  laser  speckle  imaging  [111],  diffuse  optical 
imaging  [112],  and  laminar  optical  tomography  [113])  provide  relatively  large  fields  of  view  but 
with  spatial  resolutions  of  hundreds  of  microns  to  millimeters,  microscopic  methods  (two  photon 
and  confocal  microscopy)  enable  micron-scale  resolutions.  Two-photon  microscopy  [114],  in 
particular,  is  widely  used  in  structural  and  functional  imaging  at  the  cellular  and  subcellular 
levels.  However,  two-photon  microscopy  has  imaging  speed,  penetration  depth,  and  field  of 
view. 

OCT  has  a  unique  combination  of  high  imaging  speed,  penetration  depth,  field  of  view,  and 
resolution;  and  therefore  occupies  an  important  niche  between  macroscopic  and  microscopic 
optical  imaging  technologies.  OCT  is  uniquely  capable  of  investigating  important  questions  in 
neuroscience  that  are  either  difficult  or  impossible  to  address  with  the  above  techniques.  Relative 
to  the  optical  imaging  technologies  currently  used  for  neuroscience  research,  OCT  has  several 
advantages:  1.  faster  volumetric  imaging  rates,  2.  greater  depth  penetration  with  microscopic 
resolution,  3.  quantitative  measurements  of  tissue  dynamics  (red  blood  cell  and  intra-cellular 
organelle  dynamics),  and  4.  the  ability  to  use  intrinsic  contrast,  rather  than  exogenous  fluorescent 
markers. 
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During  the  previous  funding  cycle,  we  developed  a  powerful  array  of  OCT  methods  for  in 
vivo  imaging  of  vascular  and  cellular  dynamics  in  living  animal  brains.  Specifically,  we 
developed  and  validated  Doppler  algorithms  for  absolute  CBF  measurements  in  arteries  and 
veins,  which  do  not  require  measuring  the  vessel  Doppler  angles  [18].  We  also  developed  rapid 
volumetric  imaging  of  blood  perfusion  patterns  (i.e.,  angiograms)  with  OCT  [115],  which 
enables  visualization  of  perfusion  not  only  in  arteries  and  veins,  but  also  in  capillaries.  These 
technologies  were  combined  and  applied  for  a  study  of  how  CBF  is  affected  by  ischemic  stroke 
in  animal  models  [116].  We  then  integrated  OCT  with  dynamic  light  scattering  for  microscopic 
imaging  of  tissue  dynamics  (DLS-OCT)  [117].  This  technology  enabled  us  to  simultaneously 
obtain  3D  maps  of  the  axial  and  transverse  flow  velocities  and  diffusion  coefficient,  enabling  in 
vivo  monitoring  of  not  only  blood  perfusion,  but  also  intra-cellular  organelle  motility  within 
individual  cells  [118].  Finally,  we  developed  OCT  techniques  to  quantitatively  measure  blood 
flow  within  the  capillary  network.  We  discovered  that  individual  RBC  passages  can  be  captured 
by  dynamic  OCT  imaging  and  utilized  to  measure  the  RBC  speed,  flux  and  density  over 
hundreds  of  capillaries  [119].  These  OCT  technologies  developed  for  neuro imaging  in  the 
previous  funding  cycle  lead  to  technical  advances  that  will  enable  us  to  study  a  range  of 
questions  which  could  not  be  addressed  with  existing  technologies.  This  series  of  advances  will 
also  enable  a  number  of  novel  experiments  investigating  neurovascular  coupling,  interpretation 
of  human  fMRI  data,  and  cerebrovascular  pathologies  including  those  related  to  brain  trauma. 
These  advances,  though  developed  for  neuroimaging,  can  be  translated  to  other  fields  such  as 
retinal  imaging  or  tumor  biology. 

4.2  Quantitative  Measurements  of  Absolute  Cerebral  Blood  Flow  (CBF) 

In  the  normal  brain,  cerebral  blood  flow  (CBF)  is  regulated  to  satisfy  metabolic  demand 
determined  by  brain  activity.  CBF  is  required  for  the  supply  of  oxygen  and  nutrients  and  the 
removal  of  metabolites  and  heat.  Alterations  in  CBF  may  contribute  to  the  pathogenesis  of 
stroke,  trauma,  and  Alzheimer's  disease  [120],  The  gold  standard  for  determination  of  regional 
blood  flow  is  autoradiography.  Although  autoradiographic  methods  provide  three-dimensional 
spatial  information,  they  cannot  assess  the  temporal  evolution  of  CBF  changes  [121].  Therefore, 
studies  of  disease  progression  or  response  to  treatment  using  intra-animal  comparisons  cannot  be 
performed,  limiting  the  utility  of  this  method. 

Methods  of  flow  measurement  based  on  magnetic  resonance  imaging  [122]  and  positron 
emission  tomography  [123]  provide  spatial  maps  of  CBF,  but  are  limited  in  their  temporal  and 
spatial  resolution.  Two-photon  microscopy  can  calculate  both  transverse  velocity  and  linear 
density  of  RBCs  [124],  which  may  be  combined  to  calculate  RBC  flux.  However,  the 
determination  of  CBF  requires  knowledge  of  the  flow  profile  as  well  as  vessel  orientation,  which 
can  be  difficult  to  measure.  Scanning  laser  Doppler  techniques  can  be  used  to  obtain  spatially- 
resolved  relative  CBF  images  by  moving  a  beam  across  the  field  of  interest  [125].  Laser  speckle 
imaging  [126]  has  recently  been  applied  for  high  spatio temporal  resolution  relative 
measurements  of  CBF  [111].  However,  it  is  difficult  to  relate  correlation  times  (laser  speckle)  or 
the  width  of  the  power  spectral  density  (laser  Doppler)  to  absolute  velocities  and  flow  of  red 
blood  cells,  because  the  number  of  moving  particles  that  the  light  interacts  with  and  their 
orientations  are  unknown  [127].  Doppler  OCT  overcomes  many  of  the  limitations  of  these 
aforementioned  optical  techniques  by  using  path  delay  to  gate  out  multiply  scattered  light  [104, 
128-130].  Using  “coherence  gating”  [131],  OCT  enables  measurement  of  velocity  averaged  over 
well-defined  volumes,  detennined  by  the  spot  size  in  the  transverse  direction  and  coherence 
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length  in  the  axial  direction.  Thus,  Doppler  OCT  avoids  the  limitations  associated  with  laser 
Doppler  and  laser  speckle  imaging. 

The  MGH  group  has  developed  methods  and  algorithms  for  the  quantification  of  absolute 
CBF  using  Doppler  OCT.  Figure  27  (B,  C)  show  an  en  face  image  of  the  z  projection  of  velocity 
in  the  cortex  at  the  depths  of  100  um  and  25  um,  respectively.  As  shown  in  Figure  27  (A),  by 
integrating  the  flow  profile  over  a  vessel  cross-section,  it  is  possible  to  obtain  the  flow  through 
that  vessel.  Notably,  this  method  does  not  require  the  explicit  measurement  of  the  vessel  angle 
relative  to  the  probe  beam.  Using  this  en  face  Doppler  OCT  flow  measurement,  we  compared  the 
results  of  OCT  flow  measurements  to  measurements  using  the  gold  standard  of  hydrogen 
clearance  [132].  The  results  are  shown  in  Figure  27  (C)  and  demonstrate  that  en  face  Doppler 
OCT  absolute  blood  flow  measurements  correlate  with  hydrogen  clearance  absolute  blood  flow 
measurements. 
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Figure  27.  Doppler  OCT  for  the  quantitative  study  of  cerebrovascular  physiology.  (A)  En  face  Doppler 
principle  of  calculating  flow  without  requiring  vessel  angle.  The  vessel  area  in  the  en  face  (xy)  transverse 
plane  is  inversely  proportional  to  cos(^),  while  the  z-projection  of  the  velocity  is  proportional  to  cos(^). 
Therefore,  absolute  flow  can  be  measured  by  integrating  the  axial  Doppler  velocity  over  the  vessel  area. 
The  figure  shows  an  example  of  a  venule  showing  the  z  projection  of  velocity  over  the  vessel  cross-section 
(upsampled  to  reduce  pixellation).  (B)  En  face  DOCT  image  showing  axial  projection  of  velocity  as  a 
function  of  transverse  position.  (C)  Zoomed  DOCT  images  of  boxed  region  at  a  25  pm  depth.  (D,  E) 
Validation  of  DOCT  absolute  CBF  using  simultaneous  hydrogen  clearance  in  a  rat  thinned  skull 
preparation.  (D)  OCT  angiogram  showing  the  outline  of  the  hydrogen  clearance  electrode  and  field  of  view 
for  absolute  flow  measurements  as  a  box,  where  the  dural  and  superficial  cortical  vessels  are  colored 
orange  and  the  deeper  cortical  regions  are  colored  green.  (E)  Plot  of  DOCT  absolute  flow  and  hydrogen 
clearance  flow  in  n=12  rats  (p  <  0.001). 
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4.3  DLS-OCT  Imaging  of  CBF  and  Cellular  Viability 


We  developed  an  integrated  dynamic  light  scattering  and  OCT  technique  for  high-resolution 
3D  imaging  of  heterogeneous  diffusion  and  flow  [117].  DLS  analyzes  fluctuations  in  light 
scattered  by  particles  to  measure  particle  diffusion  or  flow,  while  OCT  uses  coherence  gating  to 
collect  only  light  scattered  from  a  small  volume  for  high  resolution  structural  imaging. 
Therefore,  the  integration  of  DLS  and  OCT  enables  high  resolution  3D  imaging  of  diffusion  and 
flow.  We  derived  a  rigorous  theory  for  the  field  autocorrelation  function  of  the  complex-valued 
OCT  signal,  where  static  and  moving  particles  are  mixed  within  the  OCT  resolution  volume  and 
the  moving  particles  can  exhibit  either  diffusive  or  translational  motion.  Based  on  this  theory,  we 
developed  a  fitting  algorithm  to  estimate  dynamic  parameters  including  the  axial  and  transverse 
velocities  and  the  diffusion  coefficient  at  every  voxel  of  OCT.  We  validated  DLS-OCT 
measurements  of  diffusion  and  flow  through  numerical  simulations  and  phantom  experiments 
(Figure  28  (A)).  This  method  could  image  both  the  axial  and  transverse  velocities  of  CBF  to 
obtain  three-dimensional  maps  of  the  diffusion  coefficient  quantifying  non-translational  motions. 

In  order  to  validate  this  method,  we  performed  DLS-OCT  imaging  of  the  rodent  cerebral 
cortex,  to  measure  3D  maps  of  the  flow  velocity  and  diffusion  coefficient  (Figure  28  (B))  [118]. 
In  the  diffusion  map,  we  observed  high  diffusion  spots,  whose  locations  correspond  to  neuronal 
cell  bodies  and  whose  diffusion  coefficient  agreed  with  values  from  in  vitro  intracellular 
organelles  motion  reported  in  the  literature  (Figure  28  (C)).  This  technology  enables  label-free 
and  in  vivo  imaging  of  the  diffusion-like  motion  of  intracellular  organelles  (so  called  intracellular 
motility;  IM).  We  performed  preliminary  studies  to  demonstrate  this  method  to  monitor  CBF  and 
IM  during  a  brief  ischemic  stroke,  where  we  observed  an  induced  persistent  reduction  in  IM, 
despite  the  recovery  of  CBF  after  stroke.  As  the  preliminary  data  suggests,  DLS-OCT  provides  a 
unique  means  to  simultaneously  monitor  CBF  and  individual  cell  viability  in  vivo. 
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Figure  28.  Quantitative  imaging  of  cerebral  blood  flow  velocity  and  intracellular  motility  using 
dynamic  light  scattering  -  optical  coherence  tomography  (DLS-OCT).  (a)  Phantom  validation  of  DLS- 
OCT  measurements  of  flow  velocities  (left  and  center)  and  diffusion  coefficient  (right).  The  gray  line 
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in  the  right  panel  shows  the  Einstein-Stokes  equation.  Data  are  presented  in  mean±s.d.  (b)  DLS-OCT 
imaging  of  the  cerebral  cortex.  An  OCT  angiogram  and  DLS-OCT  maps  (absolute  velocity  and  axial 
velocity)  are  presented  as  en  face  images  using  a  maximum  projection  over  a  200  um  depth.  A 
negative  velocity  (blue)  in  the  axial  velocity  map  means  that  blood  flows  toward  the  cortical  surface. 

The  scale  bar  is  200  um.  (c)  DLS-OCT  imaging  of  neuronal  intracellular  motility.  The  velocity  and 
diffusion  images  (red  and  yellow)  are  superimposed  with  the  map  of  diffusion  with  high  coefficient  of 
determination  R2  (green).  White  circles  are  collocated  to  show  the  spatial  correlation  between  the 
positions  of  neuronal  cell  bodies  (dark  spots  in  minimum  intensity  projection,  MinIP)  and  neuronal 
intracellular  motility  (green  spots).  Not  all  cell  bodies  are  marked.  BF,  blood  flow;  IM,  intracellular 
motility;  and  VW,  vessel  wall.  Scale  bar:  200  um. 

4.2  OCT  Imaging  of  Capillary  RBC  Flow  Properties 

OCT  has  enabled  quantification  of  blood  flow  in  arteries  and  veins  [133],  and  visualization 
of  microvasculature  [11].  However,  in  order  to  quantify  blood  flow  in  the  microvasculature,  it  is 
important  to  note  that  capillaries  have  different  flow  characteristics  than  arteries  and  veins.  RBCs 
flow  one-by-one  in  capillaries,  sometimes  with  large  separations,  and  the  hematocrit  or  linear 
RBC  density  varies  significantly  across  segments  in  capillary  networks  [124].  Doppler  OCT  has 
been  widely  used  to  quantify  blood  flow  speed  in  arteries  and  veins  [133,  134],  but  its 
application  in  capillaries  is  challenging  owing  to  the  discrete  nature  of  capillary  RBC  flow  and 
the  fact  that  many  capillaries  lie  nearly  along  the  transverse  direction.  Several  decorrelation 
based  methods  have  been  proposed  to  image  the  speed  or  flow  even  in  the  transverse  direction 
[119,  135-138],  but  those  approaches  only  have  been  validated  in  large  vessels  or  glass 
capillaries  whose  diameters  are  larger  than  100  um  [136,  138-140],  Cerebral  capillaries  are  ~8 
um  in  diameter,  such  that  RBCs  only  can  pass  one-by-one.  A  second  issue  in  quantifying 
capillary  flow  properties  is  the  varying  hematocrit.  While  arteries  and  veins  have  a  relatively 
constant  hematocrit  equal  to  the  systemically  measured  hematocrit,  the  hematocrit  or  RBC  linear 
density  is  known  to  vary  significantly  across  segments  in  capillary  networks  [119,  124].  Many  of 
the  Doppler  and  decorrelation  based  methods  validated  their  measurements  of  the  RBC  speed 
using  blood  samples  with  a  constant  hematocrit,  but  the  signal  decorrelation  rate  is  generally 
hematocrit  dependent. 

Techniques  which  involve  identification  of  individual  RBC  passage  will  be  relatively  free 
from  the  above  challenges  of  the  discrete  RBC  flow  speed  and  varying  hematocrit,  and  are  able 
to  quantify  the  RBC  flux  [RBC/s].  Two-photon  microscopy  line  scanning  methods  [124]  and 
fluorescence  video  microscopy  methods  [141]  identify  individual  RBC  passage,  but  they  monitor 
only  one  or  a  few  capillaries  at  a  time.  We  developed  and  investigated  a  new  OCT  technique  that 
is  based  on  identification  of  individual  RBC  passage,  but  enables  simultaneous  measurements 
over  many  capillaries  [119].  According  to  the  Mie  scattering  theory  which  suggests  that  1  um 
wavelength  light  scattering  is  sensitive  to  scatterers  of  0.1  um  to  10  um  in  size,  high 
backscattering  (i.e.,  large  OCT  signals)  will  result  from  RBCs  as  well  as  cellular  constituents  and 
membranes.  We  validated  the  hypothesis  that  the  OCT  signal  at  a  given  voxel  is  sensitive  to  the 
transit  of  a  single  RBC  through  the  voxel  (Figure  29  (A)).  This  finding  enabled  us  to  measure 
the  RBC  speed,  flux,  and  density  over  a  large  number  of  capillaries  in  a  relatively  short  period  of 
time.  The  OCT  method  was  compared  with  the  traditional  stripe  imaging  based  technique 
(Figure  29  (B)).  This  new  OCT  method  provides  a  means  to  monitor  microvascular  flow 
dynamics  over  many  capillaries  located  at  different  depths  at  the  same  time  (Figure  29  (C)). 
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Figure  29.  Multiple-capillary  measurement  of  RBC  speed,  flux,  and  density  with  OCT.  (A)  CCD  image  of 
the  cortical  surface  through  the  cranial  window.  Scale  bar:  500  um.  (B)  The  en  face  maximum  intensity 
projection  (MIP)  of  the  3D  OCT  angiogram  over  a  0-400  um  depth.  Ten  volumes  were  averaged.  The  red 
line  indicates  the  scan  line  for  cross-sectional  dynamic  OCT  imaging.  Scale  bar:  100  um.  (C)  The  cross- 
section  of  the  3D  angiogram.  Data  are  presented  over  four  selected  capillaries  (color  circles).  Scale  bar  = 
100  um.  (D)  RBC  passage  captured  in  OCT  intensity  time  courses.  The  lines  present  the  time  courses  of 
relative  intensity  changes  at  the  four  capillaries,  with  colors  matching  to  those  in  panel  C.  The  peaks  were 
localized  over  a  few  neighboring  voxels.  Gaussian  fits  obtained  by  our  algorithm  were  overlaid  (black 
pieces).  Fitting-based  estimations  of  the  RBC  speed  (mm/s)  and  flux  (RBC/s)  are  presented  on  the  right 
side.  (E)  Comparison  between  the  fitting-  and  stripe-based  measurements  of  the  RBC  speed.  (F-H) 
Estimated  RBC  speed,  flux,  and  density  are  presented  as  color  spots  on  the  MIP  angiogram.  We  performed 
the  above  analysis  of  the  OCT  intensity  time  courses  for  each  cross-sectional  plane  in  the  volume  of 
interest.  Data  are  presented  with  relatively  comparable  ranges  (median±40%)  for  all  flow  properties. 


4.5  Rapid  Volumetric  Imaging  of  Capillary  Network  Flux 


Interest  in  blood  flow  regulation  in  the  brain  has  been  evolving  towards  understanding  the 
role  of  the  spatio-temporal  dynamics  of  capillary  networks.  In  contrast  to  arterioles,  capillaries 
have  been  reported  to  exhibit  highly  heterogeneous  responses  to  neural  activation,  capillary  by 
capillary  [142,  143],  and  nearly  stochastic  distributions  during  baseline  [141],  potentially 
masking  neural  activity-induced  responses  within  single  capillaries  [124],  Thus,  in  order  to 
characterize  capillary  network  flow  responses  to  brain  activation  with  sufficiently  high  statistical 
power,  it  is  necessary  to  measure  blood  flow  properties  over  hundreds  of  capillaries  at  the  same 
time  and  with  high  temporal  resolution  of  ~1  second  during  functional  activation.  Although  the 
technique  described  in  the  previous  section  is  faster  than  two-photon  microscopy,  it  still  requires 
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hundreds  of  seconds  to  obtain  a  map  over  a  capillary  network  consisting  of  hundreds  of 
capillaries. 

We  have  obtained  preliminary  results  using  a  novel  method  to  quantify  capillary  RBC  flux 
with  higher  temporal  resolution  using  rapid  volumetric  imaging  data,  where  only  two  scans  are 
repeated  for  each  plane.  To  achieve  this,  we  utilize  a  statistical  intensity  variation  (SIV)  metric 
quantifying  the  normalized  OCT  intensity  difference  between  the  two  consecutive  scans.  If  there 
are  no  RBCs  flowing,  the  signal  difference  is  zero.  Also,  if  there  is  no  RBC  present  during  the 
two  consecutive  scans,  the  signal  difference  is  zero.  As  the  flux  of  RBCs  increases,  the 
probability  of  an  RBC  being  present  for  one  scan,  but  not  the  other  scan  increases  and  the  signal 
difference  increases  on  average  as  the  flux  increases.  As  described  in  reference  [144]  and  show 
in  Figure  30  (A,  B),  we  have  shown  that  the  mean  SIV  is  proportional  to  the  RBC  flux  through 
both  simulations  and  measurements.  From  rapidly  scanned  volume  data,  we  used  Hessian  matrix 
analysis  to  vectorize  a  segment  path  of  each  capillary  and  estimate  its  flux  from  the  mean  of  the 
SIVs  gathered  along  the  path.  Repeating  this  process  led  to  a  3D  flux  map  of  the  capillary 
network  as  shown  in  Figure  30  (C),  which  required  only  a  10  second  data  acquisition.  Therefore, 
by  properly  splitting  the  imaging  area  and  merging  data,  this  technology  enabled  us  to  trace  the 
RBC  flux  changes  over  hundreds  of  capillaries  with  a  temporal  resolution  of  ~1  second  during 
functional  activation  (Figure  30  (D)).  This  high-throughput  monitoring  of  capillary  RBC  flux 
dynamics  will  enable  us  to  investigate  many  processes  which  were  not  possible  to  address  using 
previous  technologies.  We  believe  that  many  of  these  advanced  imaging  techniques  will  be 
translatable  to  ophthalmology  and  other  clinical  applications. 
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Figure  30.  Statistical  intensity  variation  (SIV)  analysis  for  rapid  volumetric  imaging  of  capillary  network 
flux.  (A)  Numerical  simulation  of  the  relation  between  the  mean  SIV  and  the  RBC  flux.  We  performed 
numerical  simulations  to  validate  the  mean  SIV  as  a  function  of  the  RBC  flux.  First,  we  numerically 
synthesized  OCT  time  courses  embedding  RBC  passage  peaks,  for  various  flow  speeds  and  mean  RBC 
separations.  From  each  time  course,  we  sampled  a  set  of  SIV  values  with  the  time  gap  equal  to  the  one  that 
was  used  in  experiments  (11  ms).  Then,  we  averaged  the  set  of  SIV  values  into  a  mean  value  and  compared 
it  to  an  RBC  flux  obtained  by  counting  the  RBC  passage  peaks  during  the  time  course.  The  mean  SIV  was 
proportional  to  the  RBC  flux  regardless  of  various  speeds  and  densities  as  hypothesized  (r  =0.97).  (B)  We 
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also  experimentally  investigated  the  relation  of  SIV  to  RBC  flux.  We  obtained  OCT  intensity  time  courses 
of  capillary  centers  (4  s  in  length,  n  =  22)  that  exhibited  individual  RBC  passage  as  in  the  publication  [119]. 
We  divided  each  time  course  into  four  1-s  segments;  for  each  segment,  we  counted  the  number  of  peaks  for 
true  RBC  flux  values  and  obtained  the  mean  of  randomly  sampled  SIVs;  and  then  presented  their  mean  and 
standard  deviation  (SD)  as  a  single  data  point  in  this  figure.  We  experimentally  found  that  the  mean  SIV 
was  linearly  related  with  the  RBC  flux  (r  =  0.78).  The  slope  and  intersect  will  depend  on  the  measurement 
system  and  scanning  protocol  including  the  time  gap  (At).  The  slope  was  slightly  larger  in  the  experiment, 
likely  because  the  experimental  data  had  non-zero  fluctuations  even  when  no  RBC  passes.  The  empirical 
relation  with  our  system  and  At  =  11  ms  was  (SIV)  =  0.012  +  0.032  x  (Flux).  Since  isotropic  voxels  were 
used,  this  relation  is  in  principle  independent  of  the  direction  of  RBC  flow  or  the  angle  between  the 
capillary  axis  and  the  scanning  plane.  (C)  3D  flux  map  of  the  capillary  network.  Using  our  algorithms  to 
vectorize  capillary  segment  paths  and  the  SIV  volume  data,  the  RBC  flux  of  each  capillary  segment  was 
obtained  from  the  mean  SIV  averaged  along  its  path.  This  technique  enabled  us  to  measure  the  RBC  flux 
over  hundreds  of  capillaries,  leading  to  a  3D  flux  map  of  the  capillary  network.  The  bar  graph  in  the  right 
shows  the  RBC  flux  histogram  of  178  capillaries.  (D)  Time  courses  of  RBC  flux  changes  196  capillaries 
during  functional  activation.  We  performed  dynamic  SIV  imaging  over  an  ROI  near  the  center  of 
functional  activation  in  the  somatosensory  cortex  corresponding  to  forepaw  stimulation.  A  temporal  series 
of  20  SIV  volumes  was  obtained  by  averaging  over  ten  runs  of  stimulation,  leading  to  a  temporal  series  of 
20  capillary  network  flux  maps.  It  was  possible  to  trace  how  the  RBC  flux  varies  with  activation  over  196 
capillaries.  A  change  in  the  mean  flux  averaged  across  capillaries  is  presented  by  the  thick  black  curve.  The 
peak  change  was  2.2%  and  highly  significant  (p  <  10  s).  The  black  bar  in  the  bottom  indicates  the  duration 
of  forepaw  stimulation  (3  Hz  for  T=  0  -  4  s). 


5.0  Endoscopic  OCT  Technology  and  Applications 

Barrett’s  esophagus  (BE)  is  a  common  condition  associated  with  gastroesophageal  reflux 
disease  (GERD)  in  which  the  nonnal  squamous  epithelium  of  the  esophagus  is  replaced  by 
metaplastic  columnar  epithelium  [145,  146].  In  the  US,  about  5.6%  of  the  population  is  affected 
by  this  disease  [147].  Patients  with  BE  have  an  increased  risk  of  developing  esophageal 
adenocarcinoma  [148],  which  is  a  lethal  disease  with  less  than  20%  five-year  survival  rate  [149]. 
The  incidence  of  esophageal  adenocarcinoma  (EAC)  has  increased  dramatically  over  the  past 
decades  [150,  151]. 

For  over  15  years  we  have  been  working  in  collaboration  with  Dr.  Hiroshi  Mashimo,  MD, 
PhD,  at  the  VA  Boston  Healthcare  System  and  Harvard  Medical  School  to  develop  and  apply 
advanced  endoscopic  OCT  technology  to  improve  the  management  of  BE  and  related  diseases  of 
the  gastrointestinal  (GI)  tract.  During  the  previous  grant  period  we  made  significant  advances 
towards  developing  endoscopic  OCT  as  a  clinically  accepted  imaging  modality.  We  performed 
clinical  studies  on  the  natural  history  of  BE  disease,  as  well  as  identifying  markers  that  can  be 
useful  in  planning  endoscopic  interventions  or  predict  treatment  outcomes.  In  another  set  of 
studies,  we  continued  advancing  the  state  of  the  art  of  the  endoscopic  OCT  technology  by 
developing  OCT  imaging  engines  with  ultrahigh  imaging  speeds  and  OCT  imaging  catheters 
with  high  speed,  precision  beam  scanning.  This  new  generation  ultrahigh  speed  endoscopic  OCT 
technology  enables  three  dimensional,  depth  resolved  visualization  mucosal  surface  patterns  and 
microvasculature,  promising  to  enhance  the  sensitivity  and  specificity  of  OCT  for  identifying 
early  neoplastic  changes  as  well  as  assessing  treatment  response. 

5.1  Clinical  Studies  with  Endoscopic  OCT 

Radiofrequency  ablation  (RFA)  is  an  emerging  endoscopic  therapy  for  treating  Barrett’s 
esophagus  (BE)  [40-44],  Utilizing  electrode  arrays,  RFA  catheters  deliver  radiofrequency  energy 
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to  ablate  BE  with  a  low  risk  of  complications  such  as  stricture  [45,  46].  Recent  studies  have 
shown  that  complete  eradication  of  intestinal  metaplasia  (CE-IM)  was  achieved  in  93%  of 
patients  with  dysplasia  after  two  years  [152]  and  in  92%  of  patients  with  non-dysplastic  BE 
(NDBE)  at  5-year  follow-up  [153].  However,  recurrence  of  intestinal  metaplasia  (IM)  has  been 
observed  in  13%  [152]  to  26%  [154]  of  patients  at  1  year  after  CE-IM.  Repeated  RFA  treatments 
are  generally  required  to  achieve  CE-IM  [47-49].  On  average,  CE-IM  was  achieved  after  3.4 
RFA  sessions  for  patients  with  NDBE  [153]  and  over  3.5  sessions  for  patients  with  dysplasia 
[43,  152].  After  RFA  treatment,  new  squamous  epithelium  generally  replaces  the  BE  epithelium. 
However,  residual  BE  glands  can  become  "buried"  by  the  neosquamous  epithelium  and  are  not 
visible  by  conventional  endoscopy,  even  using  chromogen  or  narrow  band  imaging  (NBI). 
Because  areas  of  neosquamous  epithelium  are  not  routinely  biopsied  in  clinical  practice,  these 
buried  glands  are  likely  under-appreciated  with  the  current  surveillance  protocol  [43,  153,  155]. 

OCT  has  the  capability  to  visualize  subsurface  features  in  the  human  esophagus  with  1-2  mm 
imaging  depth,  enabling  evaluation  of  morphology  under  the  squamous  epithelium.  The 
development  of  high  speed  endoscopic  OCT  technology  [156-159],  enables  three-dimensional 
(3D)  comprehensive  imaging  of  the  epithelium.  Because  buried  glands  cannot  be  detected  using 
standard  endoscopic  imaging,  OCT  uniquely  enables  depth-resolved  imaging  of  a  broad  area 
with  near-microscopic  resolution  and  is  a  promising  technology  for  identifying  and 
characterizing  buried  glands  before  and  after  ablative  therapies.  In  addition,  depth-resolved 
imaging  can  also  characterize  the  thickness  of  different  mucosal  layers  in  the  BE  region,  which 
may  guide  treatment  dosage  and  improve  efficacy. 

During  the  previous  contract  period,  we  performed  several  clinical  studies  of  RFA  using  a 
prototype  3D-OCT  imaging  system  developed  by  LightLab  Imaging  (a  subsidiary  of  Saint  Jude 
Medical)  combined  with  custom  built  imaging  probes.  The  SS-OCT  system  used  a  Fourier 
domain  mode-locked  laser  operating  at  60  kHz  A-scan  sweep  rate  and  1310  nm  wavelength.  The 
system  had  ~5  um  axial  and  ~15  um  transverse  resolution  in  tissue.  The  OCT  imaging  probe  had 
a  2.5  mm  outer  diameter  and  was  radially  scanned  at  60  Hz  and  pulled  back  at  1  mm/s.  3D-OCT 
data  sets,  covering  a  region  of  8  mm  (circumference)  x  20  mm  (length)  x  2  mm  (depth)  were 
acquired  in  20  seconds. 

In  one  study,  we  quantified  buried  glands  using  endoscopic  3D-OCT  in  a  cohort  of  27 
patients  undergoing  RFA  treatment  [160].  In  a  second  study,  we  identified  structural  markers 
that  predicted  RFA  treatment  response  using  endoscopic  3D-OCT  in  a  cohort  of  33  patients 
undergoing  RFA  treatment  [50].  Results  of  these  two  studies  will  be  summarized  in  the  next  two 
sections. 

Characterization  of  Buried  Glands  Pre-  and  Post-  Radio  frequency  Ablation  using  3D-OCT 

Figure  31  shows  examples  of  cross-sectional  OCT  images  (A,  B)  and  corresponding 
histology  (c)  of  the  gastroesophageal  junction  (GEJ).  Buried  glands  (red  arrows)  were  identified 
with  3D-OCT  at  the  GEJ  of  a  patient  during  the  imaging  procedure.  Biopsies  were  performed  on 
landmarks  identified  during  OCT  imaging  using  video  endoscopy.  The  histological  micrograph 
shown  in  Figure  31  (C)  confirmed  the  presence  of  buried  glands  with  intestinal  metaplasia 
underneath  neosquamous  epithelium. 
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Figure  30.  Cross-sectional  OCT  images  and  corresponding  histology  showing  buried  glands  (red  arrows) 
from  one  patient  in  vivo  (a-c).  Histological  micrograph  shown  in  (c)  confirmed  the  presence  of  buried 
glands  with  intestinal  metaplasia.  SE,  squamous  epithelium.  Scale  bars:  500  um. 

From  the  entire  patient  cohort,  buried  glands  were  identified  in  72%  (13/18)  of  patients  in  the 
pre-CE-IM  group  and  63%  (10/16)  of  patients  in  the  post-CE-IM  group.  Overall,  620  and  114 
buried  glands  were  marked  in  3D-OCT  datasets  from  patients  in  the  pre-CE-IM  and  post-CE-IM 
groups,  respectively.  For  comparison,  Figure  32  shows  histograms  of  the  gland  sizes  (Figure  32 
(A))  and  the  distance  from  the  squamous-columnar  junction  (SCJ)  (Figure  32  (B)).  On  average, 
the  size  of  the  buried  glands  was  214±126  um  in  the  pre-CE-IM  group,  and  219±1 1 1  um  in  the 
post-CE-IM  group,  respectively.  Furthermore,  the  distribution  of  the  buried  glands  relative  to  the 
SCJ  was  also  similar  in  both  groups  (Figure  32  (B)).  67%  (415/620)  and  69%  (79/114)  of  the 
glands  were  observed  within  5  mm  from  the  SCJ  in  the  pre-CE-IM  group  and  post-CE-IM  group, 
respectively.  The  average  distance  of  these  glands  from  the  SCJ  was  3.6±2.7  mm  in  the  pre-CE- 
IM  group  and  3.8±3.9  mm  in  the  post-CE-IM  group  and  was  not  statistically  significant  between 
the  two  groups  (p  =  0.54). 
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Figure  32.  Histograms  of  (a)  the  buried  gland  size  and  (b)  distance  from  the  SCJ  in  the  Pre-CE-IM  and 
Post-CE-IM  groups.  No  significant  difference  was  observed  in  either  the  gland  size  or  gland  distribution 
between  the  two  groups. 

This  study  demonstrates  that  3D-OCT  provides  sufficient  sensitivity  to  detect  buried  glands 
pre-  and  post-CE-IM.  Interestingly,  a  very  high  prevalence  of  buried  glands  was  observed  in  both 
groups,  with  72%  of  patients  pre-CE-IM  and  63%  in  the  post-CE-IM  group.  The  number  of 
buried  glands  per  patient  decreased  significantly  after  CE-IM  is  achieved  (p  =  0.02).  However, 
the  size  and  distribution  of  the  buried  glands  relative  to  the  SCJ  in  patients  with  short  segment 
BE  did  not  change  significantly  post-CE-IM. 
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Endoscopic  3D-OCT  Markers  of  Radio  frequency  Ablation  Treatment  Response 

Figure  33  (A)  shows  an  example  endoscopic  image  of  the  gastroesophageal  junction  (GEJ) 
before  focal  RFA  treatment.  3D-OCT  imaging  was  performed  at  the  GEJ  over  the  area  of  BE 
(pink  colored).  Figure  33  (B)  shows  a  representative  cross-sectional  OCT  image  obtained  at  the 
GEJ  prior  to  the  RFA  treatment,  together  with  histology  obtained  from  the  same  site  (Figure  33 
(C)).  The  thickness  of  the  BE  epithelium  was  measured  from  the  OCT  image  as  distance 
between  the  BE  surface  to  the  top  of  the  lamina  propria  (LP)  /  muscularis  mucosa  (MM)  layer. 


Figure  33.  (a)  Representative  endoscopic  image  of  the  GEJ  before  RFA  treatment,  (b)  Representative 
cross-sectional  endoscopic  OCT  image  and  (c)  corresponding  histology  illustrating  the  BE  epithelium 
thickness  measurement. 

The  BE  thickness  measured  with  OCT  prior  to  RFA  was  found  to  be  a  predictor  of  the  RFA 
treatment  response.  Both  the  average  and  maximum  BE  thickness  prior  to  RFA  were 
significantly  thinner  for  the  CE-IM  group  compared  to  the  non-CE-IM  group  (average  BE 
thickness,  257±60  um  versus  403±86  um,  p<0.0001).  The  scatter  plot  in  Figure  34  (A)  shows 
the  difference  between  the  average  BE  thickness  for  the  two  groups.  Figure  34  (B)  shows  ROC 
curves  using  average  and  maximum  BE  thickness  to  predict  RFA  treatment  response.  The  area- 
under-the-curve  (AUC)  was  0.942  (p  <  0.001)  and  0.934  (p  <  0.001)  using  the  average  and 
maximum  BE  thickness,  respectively.  An  average  BE  thickness  of  333  um  was  determined  from 
the  ROC  curve  to  achieve  the  best  prediction  accuracy.  Using  this  decision  threshold,  a 
sensitivity  of  92.3%  (12/13),  specificity  of  85%  (17/20),  PPV  of  80%,  NPV  of  94.4%  and  an 
accuracy  of  87.9%  (29/33)  was  obtained  for  predicting  treatment  response  evaluated  by  the 
presence  or  absence  of  endoscopically  visible  residual  BE  at  follow-up  visit  using  the  average 
BE  thickness  measured  with  OCT  prior  to  the  RFA  treatment.  The  BE  thickness  measured  with 
OCT  was  not  correlated  with  the  length  of  the  BE  (p  =  0.88)  or  the  number  of  prior  RFA 
treatments  (p  =  0.24). 
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Figure  34.  (A)  Scatter  plot  of  the  average  BE  epithelium  thickness  measured  by  OCT.  Blue  circles:  Non- 
CE-IM  group;  Red  crosses:  CE-IM  group;  Green  dotted  line:  discrimination  threshold  at  333  um  as 
determined  from  the  average  BE  thickness  ROC  curve  in  (B).  (B)  ROC  curves  of  treatment  response 
prediction  using  average  (green)  and  maximum  (blue)  BE  thickness.  The  area-under-the-curve  (AUC) 
values  were  0.942  (p  <  0.001)  and  0.934  (p  <  0.001)  using  the  average  and  maximum  BE  thickness, 
respectively. 

In  summary,  these  studies  show  that  endoscopic  3D-OCT  can  identify  structural  markers, 
such  as  BE  thickness  and  residual  glandular  structure  which  are  predictors  of  RFA  treatment 
response  at  follow  up.  3D-OCT  may  provide  information  which  will  enable  real-time  treatment 
decisions,  improving  the  efficacy  of  RFA  treatment  and  reducing  the  number  of  required 
treatments  in  the  future. 

5.2  Endoscopic  OCT  Technology  Development 

Ultrahigh  Speed  Endoscopic  Swept  Source  OCT  System 

Our  ongoing  collaboration  with  Praevium  Research  and  Thorlabs  Inc.  in  VCSEL  light  source 
development  has  enabled  us  to  design  and  develop  ultrahigh  speed  endoscopic  swept  source 
OCT  (SS-OCT)  technology  that  is  10-20  times  faster  than  commercially  available  endoscopic 
OCT  technology.  We  have  translated  this  ultrahigh  speed  technology  to  the  endoscopy  clinic, 
working  in  collaboration  with  Dr.  Hiroshi  Mashimo,  M.D.,  Ph.D.  at  the  Boston  Veterans  Affairs 
Healthcare  System  and  Harvard  Medical  School.  Ongoing  clinical  studies  are  investigating  the 
diagnostic  potential  of  this  advanced  endoscopic  OCT  technology. 

The  ultrahigh  speed  endoscopic  SS-OCT  system  operates  at  1300  nm  using  a  MEMS  VCSEL 
light  source.  The  system  schematic  is  shown  in  Figure  35.  The  laser  sweep  is  triggered  by  a 
function  generator  at  a  frequency  of  300  kHz.  Both  forward  and  backward  sweeps  are  used  for 
imaging,  resulting  in  an  effective  A-scan  rate  of  600  kHz.  The  system  uses  a  dual  circulator 
interferometer  design  to  maximize  detection  efficiency.  A/D  sampling  is  performed  using  optical 
clocking  from  a  Mach-Zehnder  interferometer  to  sample  the  OCT  fringe  data  linear  in 
wavenumber,  greatly  simplifying  processing  and  reducing  data  volume.  Acquisition  was 
performed  using  a  high  speed  A/D  card  (AlazarTech)  and  custom  C++  software.  Axial  resolution 
was  11  um  in  air  and  8  um  in  tissue,  with  an  imaging  range  of  3.3  mm  in  air.  The  ultrahigh 
imaging  speeds  achieved  by  this  technology  enable  both  3D  volumetric  OCT  as  well  as  OCT 
angiography. 
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Figure  35.  Layout  of  ultrahigh  speed  VCSEL  swept  source  OCT  system  for  endoscopic  applications.  C: 
circulator.  HVA:  high  voltage  amplifier.  AWG:  arbitrary  waveform  generator.  MZI:  Mach-Zehnder 
interferometer.  P:  photodetector.  OC:  optical  clock.  DA:  differential  amplifier.  RM:  reference  mirror. 
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Figure  36.  Schematic  of  micromotor  imaging  probe.  The  distal  micromotor  provides  high  speed  and  more 
stable  beam  scanning  than  previous  proximally  actuated  torque  cable  designs.  The  probe  is  small  enough  to 
pass  through  the  3.7  mm  accessory  port  of  a  standard  endoscope. 


Micromotor  Endoscopic  Imaging  Probes 

To  fully  utilize  ultrahigh  speed  imaging,  we  have  concurrently  developed  imaging  probe 
devices  capable  of  high  speed  beam  scanning.  Using  a  micromotor  as  an  optical  beam  scanner  it 
is  possible  to  develop  an  imaging  probe  that  is  only  3.3  mm  outer  diameter  and  less  than  20  mm 
rigid  length  (Figure  36).  The  small  probe  size  is  important  because  it  can  be  inserted  into  the  3.7 
mm  accessory  channel  of  an  endoscope  and  used  as  an  adjunct  to  standard  endoscopy.  The 
micromotor  scans  in  the  rotary  direction,  and  the  probe  is  pulled  back  at  the  proximal  end  with  a 
motorized  stage  at  a  calibrated  speed  to  generate  a  3D  OCT  volume.  High  speed  distal  scanning 
avoids  the  need  for  proximal  torque  cable  rotation,  which  is  a  major  source  of  non-uniform 
rotational  distortion  in  previous  OCT  imaging  probes.  This  technology  was  validated  in  an  in 
vivo  rabbit  imaging  study  [6].  The  imaging  speed  was  1  MHz  A-scan  rate  and  400  frames  per 
second.  The  focused  transverse  spot  size  of  the  probe  output  beam  was  8  um  (FWHM). 
Volumetric  images  were  acquired  in  the  rabbit  esophagus  and  rectum,  as  well  as  ex  vivo  human 
colon  tissue.  This  study  demonstrated  that  ultrahigh  speed  OCT  and  high  speed  scanning  probe 
technology  could  acquire  image  volumes  over  large  fields  at  microscopic  resolution. 

The  instrument  was  then  engineered  for  use  in  the  endoscopy  suite  [161].  3D-OCT  imaging 
was  performed  in  patients  with  a  cross  section  of  pathologies  undergoing  upper  and  lower 
endoscopy.  The  use  of  distally  actuated  imaging  micromotor  catheters  enabled  improved  access, 
such  as  volumetric  imaging  in  the  small  intestine  and  the  assessment  of  hiatal  hernia  using 
retroflex  imaging.  The  high  rotational  scanning  stability  of  the  micromotor  enabled  3D 
volumetric  imaging  with  micron  scale  volumetric  accuracy  for  both  en  face  OCT  and  cross- 
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sectional  imaging,  as  well  as  OCT  angiography  (OCTA)  for  3D  visualization  of  subsurface 
microvasculature. 


Figure  37.  (A,  B,  C)  Volumetric  image  set  of  rabbit  esophagus  showing  en  face,  cross-sectional  and 
longitudinal  planes  respectively.  The  layered  architecture  of  the  esophageal  wall  is  clearly  resolved.  (D,  E) 
Volumetric  image  of  ex  vivo  human  colon  showing  cross-sectional  and  en  face  planes  respectively.  Surface 
architecture  of  colonic  crypts  are  visible. 


5.3  Novel  Diagnostics  with  Endoscopic  OCT 


Progression  to  esophageal  adenocarcinoma  (EAC)  is  characterized  by  a  series  of  cytological 
and  tissue  architectural  changes  that  occur  in  the  squamous  epithelium  (SE).  Following  repeated 
acid  exposure  due  to  gastroesophageal  reflux  disease  (GERD),  SE  is  initially  replaced  by 
columnar  lined  esophagus,  a  condition  known  as  Barrett’s  esophagus  (BE).  BE  in  turn  progress 
to  low  grade  (LGD)  and  later  to  high  grade  dysplasia  (HGD),  and  eventually  becomes  EAC  [39]. 
Dysplasia  of  any  grade  increases  cancer  risk  by  an  order  of  magnitude  [162],  and  HGD  is 
associated  with  a  10-60%  risk  of  developing  adenocarcinoma  within  3-5  years  [163-165]. 
Although  EAC  is  highly  lethal  once  the  cancer  arises,  prognosis  is  very  favorable  if  it  can  be 
detected  and  treated  at  a  dysplastic  state.  Treatment  methods  for  BE  with  HGD,  such  as 
radiofrequency  ablation  (RFA)  have  been  shown  to  have  high  success  rates.  Therefore,  detection 
of  dysplasia  is  of  paramount  importance  in  reducing  progression  to  EAC. 
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Over  the  past  decade  several  studies  have  assessed  the  utility  of  OCT  to  identify  dysplasia 
[166,  167].  However  limitations  in  imaging  speed  made  it  difficult  to  acquire  high  density 
volumetric  data  sets  which  are  required  in  order  to  generate  en  face  OCT  images.  Furthermore, 
previous  OCT  imaging  probes  employed  proximal  scanning  that  necessitated  the  rotation  of  a 
long  (~2  m)  torque  cable.  This  caused  non-uniform  rotational  distortion  (NURD)  from  friction 
between  the  torque  coil  and  probe  sheath,  distorting  en  face  image  features  [168].  Therefore,  to 
date,  diagnostic  criteria  used  in  clinical  OCT  studies  have  been  based  on  cross-sectional  image 
features  [166,  167].  However,  these  studies  have  shown  a  limited  diagnostic  accuracy  and 
moderate  interobserver  agreement.  The  results  of  previous  endoscopic  OCT  studies  suggest  that 
cross-sectional  OCT  alone  might  not  be  sufficiently  accurate  for  detection  of  dysplasia. 

We  have  recently  designed  an  ultrahigh  speed  endoscopic  swept  source  OCT  system  and 
beginning  in  2013  deployed  it  to  the  endoscopy  clinic  at  the  Boston  VA  Healthcare  System  and 
Harvard  Medical  School  to  perfonn  clinical  endoscopic  studies.  This  advanced  technology 
enables  new  imaging  techniques  which  promise  to  improve  the  diagnostic  potential  of 
endoscopic  OCT.  We  have  demonstrated  that  the  high  speed  imaging  and  precise  beam  scanning 
allows  depth  resolve  enface  OCT  and  visualization  of  mucosal  surface  patterns,  which  may  be 
more  sensitive  markers  of  dysplasia  than  cross-sectional  image  features.  In  addition,  we  have 
also  perfonned  the  first  demonstration  of  endoscopic  OCT  angiography  (OCTA),  enabling  three 
dimensional  visualization  of  the  microvasculature  without  the  need  of  exogenous  contrast  agents. 
This  result  was  featured  on  cover  of  Gastroenterology,  147,  6,  December  2014,  Tsai,  et  al.  This 
result  is  significant  not  only  because  it  represents  a  highly  challenging  technical 
accomplishment,  but  also  because  the  combination  of  microvascular  features  with  structural 
features  promises  to  improve  the  sensitivity  and  specificity  for  detecting  early  neoplastic 
changes.  The  next  sections  will  summarize  the  clinical  results  we  have  obtained  with  advanced 
endoscopic  OCT. 

En  Face  OCT  and  Pit  Patterns 

Conventional  endoscopy  produces  views  that  visualize  surface  morphology.  Therefore,  a 
number  of  imaging  techniques  have  been  developed  to  not  only  characterize  features  observed 
on  the  endoscopic  views  and  correlate  them  with  pathology,  but  also  to  improve  the  contrast  and 
visibility  of  the  surface  structures.  Imaging  techniques  such  as  magnifying  endoscopy, 
chromoendoscopy,  narrowband  imaging  (NBI)  and  confocal  laser  endomicroscopy  (CLE)  have 
been  extensively  studied  to  enhance  the  visualization  of  the  surface  patterns  and  improve  the 
diagnostic  yield  of  endoscopy.  However,  these  techniques  suffer  from  limitations  including  low 
diagnostic  accuracy,  high  interobserver  variation,  limited  coverage,  extended  procedural  time 
and  the  need  for  exogenous  contrast  agents;  which  so  far  hindered  their  clinical  acceptance. 

Using  ultrahigh  speed  endoscopic  swept  source  OCT  technology,  we  have  imaged  more  than 
40  patients  to  date.  Continuting  enrollment  is  ongoing.  The  ultrahigh  speed  endoscopic  OCT 
system,  combined  with  high  speed  and  precision  distal  scanning  probes  enabled  depth  resolved 
visualization  of  mucosal  surface  patterns  with  high  resolution  and  contrast.  Figure  38  shows 
representative  en  face  OCT  images  from  a  variety  of  normal  mucosal  surface  patterns 
pertinent  to  both  the  upper  and  lower  GI  tracts,  as  well  as  an  example  of  non-dysplastic 
Barrett’s  esophagus  (NDBE).  Figure  38  (A)  shows  the  colonic  folds  and  regular  circular  pit 
pattern  associated  with  nonnal  colon.  Figure  38  (B)  shows  nonnal  squamous  epithelium 
(SE)  which  shows  homogeneous  tissue  structure  as  expected.  Figure  38  (C)  is  from  an 
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area  with  NDBE  and  shows  regularly  organized  BE  pit  patterns.  Figure  38  (D)  shows 
regular  gastric  pit  patterns  associated  with  normal  stomach. 


Esophagus  (NDBE).  (A)  Colonic  mucosa.  (B)  Esophageal  mucosa.  (C)  Barrett’s  mucosa.  (D)  Gastric 


mucosa. 


Figure  39.  Representative  en  face  OCT  images  of  images  of  pathologies  in  the  upper  GI  tract.  (A)  Low 
grade  dysplasia  (LGD).  (B)  Esophageal  adenocarcinoma  (EAC),  grade  Tl.  (B,  D)  Corresponding  histology 
from  (A)  and  (C),  respectively.  Arrows  point  to  irregular  pit  pattern  architecture. 

Figure  39  shows  representative  en  face  OCT  images  of  upper  GI  pathologies  and 
corresponding  histology.  Figure  39  (A)  is  from  a  region  with  low  grade  dysplasia  (LGD),  while 
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Figure  39  (C)  shows  intramucosal  adenocarcinoma  (IMC).  Compared  with  regularly  shaped  and 
distributed  pit  patterns  of  NDBE,  LGD  and  IMC  exhibit  irregularity  and  distortion  on  the  pit 
patterns.  The  appearance  of  the  mucosal  surface  patterns  resembles  to  “ridge/villous”  and 
“irregular/distorted”  patterns  observed  in  the  endoscopy  imaging  literature,  which  are  associated 
with  the  presence  of  malignancy  [169]. 

As  demonstrated  in  Figures  40  and  41,  the  ultrahigh  speed  endoscopic  OCT  system  is 
capable  of  visualizing  mucosal  surface  pattern  with  high  contrast  and  resolution.  En  face  OCT 
images  show  marked  differences  between  non-dysplastic  and  dysplastic  BE  tissue.  We  are 
currently  continuing  patient  recruitment  as  well  as  analyzing  the  OCT  images  to  establish  a 
comprehensive  diagnostic  criteria  and  performing  a  sensitivity/specificity  test  to  validate  the 
utility  of  enface  OCT  in  detecting  early  malignancy  in  the  GI  tract. 

Endoscopic  OCTA  and  Microvasculature 

An  important  diagnostic  criteria  used  in  CLE  and  NBI  for  identification  of  malignancy  is  the 
regularity  and  distribution  of  tissue  vasculature.  However,  conventional  structural  OCT  imaging 
does  not  have  vascular  contrast.  Only  larger  vessels  can  be  visualized  as  hyposcattering 
structures  as  they  cast  shadows  to  the  underlying  tissues.  To  overcome  this  limitation,  a  set  of 
phase  and  intensity-sensitive  signal  processing  techniques  have  been  developed  that  can  produce 
vascular  contrast  from  OCT  images.  These  techniques  have  been  broadly  named  as  Doppler 
OCT  [104,  128]  and  OCTA  [11,  76,  170].  Although  OCTA  has  been  demonstrated  for 
visualizing  larger  vasculature  in  an  animal  model  [158],  speed,  resolution  and  stability  of 
previous  systems  have  not  been  sufficient  for  endoscopic  clinical  applications.  Consequently, 
angiography  imaging  studies  in  human  subjects  have  been  limited  to  single  cross-sectional  2D 
frames  [171],  making  assessment  of  GI  pathologies  challenging. 

The  ultrahigh  speed  imaging  and  stable  beam  scanning  achieved  by  ultrahigh  speed 
endoscopic  OCT  technology  makes  endoscopic  OCTA  feasible.  We  have  recently  demonstrated 
endoscopic  OCTA  from  both  upper  and  lower  GI  tracts  of  patients  undergoing  surveillance 
endoscopy  and  colonoscopy  [2].  Figure  40  shows  example  enface  OCT  and  OCTA  images  from 
nonnal  esophagus  and  recto-anal  junction.  The  structural  en  face  OCT  images  in  Figures  40  (A, 
B,  E)  show  details  of  tissue  architecture  and  mucosal  morphology,  however,  vascular 
information  is  limited  to  large  vessels.  The  corresponding  OCTA  images  shown  in  Figures  40 
(C,  D,  F),  on  the  other  hand,  visualize  microvasculature  with  high  resolution  and  contrast. 
OCTA  is  advantageous  over  CLE  since  it  does  not  require  intravenous  injection  of  exogenous 
contrast  agents  (such  as  fluorescein).  Furthermore,  it  also  allows  depth  resolved  visualization  of 
microvasculature  over  a  large  field  of  view  (16  mm  x  10  mm),  which  is  not  attainable  by  CLE 
and  NBI.  Therefore,  endoscopic  OCTA  promises  to  be  a  valuable  adjunct  to  OCT  structural 
imaging. 
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Figure  40.  Endoscopic  optical  coherence  tomography  (OCT)  and  angiography  (OCTA)  of  normal 
esophagus  (A-D)  and  recto-anal  junction  (E,  F).  (A)  En  face  OCT  image  at  240  pm  depth  (lamina  propria 
[LP]  layer).  (B)  Cross-sectional  OCT  image  along  the  pullback  direction.  (C)  En  face  OCT  angiogram  in 
the  LP  layer.  (D)  Cross-sectional  OCTA  along  the  pullback  direction.  (E)  En  face  OCT  image  at  220  pm 
depth  (LP  layer).  (F)  OCT  angiogram  shows  the  surface  vasculature  of  the  crypts  in  the  rectum.  Red  arrows, 
anal  glands;  blue  arrows,  rectal  glands;  white  bar,  1  mm. 


Figure  41.  Demonstration  of  the  NURD  correction  algorithm  for  en  face  OCT  and  OCTA  images.  En  face 
OCT  images  before  (A)  and  after  (B)  applying  the  correction  algorithm,  respectively.  OCT  angiograms 
before  (C)  and  after  (D)  applying  the  correction  algorithm,  respectively.  All  figures  are  generated  by 
projection  over  140  pm  of  depth,  and  about  150-500  pm  beneath  tissue  surface.  (A)  and  (C)  are  taken  from 
the  rectum,  and  (B)  and  (D)  are  taken  from  esophagus  of  two  patients  who  were  undergoing  endoscopic 
surveillance.  Scale  bars  are  200  pm  for  (A)  and  (C),  and  500  pm  for  (B)  and  (D). 
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Both  en  face  OCT  and  OCTA  can  have  imaging  artifacts  due  to  mechanical  and  patient- 
related  motion.  Image  processing  to  generate  OCTA  involves  calculating  the  correlation  between 
subsequently  acquired  OCT  frames  in  order  to  detect  motion  from  blood  flow;  hence  the  quality 
of  the  OCT  angiograms  is  particularly  sensitive  to  parasitic  motion  or  beam  scanning  instability. 
Motion  artifacts  include  cardiac  and  respiratory  motion,  as  well  as  non-unifonn  rotational 
distortion  (NURD)  of  the  imaging  catheter.  NURD  is  a  major  cause  of  motion  artifacts 
encountered  in  endoscopic  OCT  systems  [168].  While  using  a  micromotor  imaging  probe 
significantly  reduces  the  severity  of  NURD,  residual  non-uniformity  of  the  motor  rotation  on  the 
order  of  several  milli-radians  can  still  cause  distortion  in  the  en  face  OCT  and  OCTA  images  as 
shown  in  Figure  40  (A,  B).  To  address  this  problem,  we  have  developed  a  NURD  correction 
algorithm  that  resamples  cross-sectional  OCT  images  based  on  fiducial  markers  located  on  the 
OCT  probe  [172].  We  have  shown  that  this  method  achieves  a  ~15  fold  improvement  in 
rotational  scanning  stability  and  effectively  suppresses  artifacts  caused  by  NURD  (Figure  41  (C, 
D)). 

In  summary,  the  development  of  advanced  ultrahigh  speed  endoscopic  OCT  technology 
based  on  the  VCSEL  light  source  allowed  visualization  of  mucosal  surface  patterns  and  the 
development  of  endoscopic  OCTA  with  micron-scale  resolution  and  high  contrast.  Our 
preliminary  studies  suggest  that  these  new  capabilities  are  promising  for  improving  the  sensitivy 
and  specificity  to  detect  pathology  in  the  GI  tract  including  dysplasia  in  Barrett’s  esophagus.  The 
detection  of  dysplasia  in  Barrett’s  esophagus  has  been  an  ongoing  challenge  for  the  endoscopic 
imaging  community  and  larger  scale  prospective  studies  using  integrated  structural  OCT  and 
OCTA  are  planned  to  investigate  diagnostic  perfonnance. 

6.0  Pathology  Imaging 

Cancer  is  the  second  leading  cause  of  death  in  the  United  States,  responsible  for  577,000 
deaths  in  2012  [173].  Following  the  detection  of  malignancy  through  imaging,  routine  screening 
or  other  methods,  the  standard  of  care  for  many  types  of  cancer,  including  malignant  tumors  of 
the  breast  [174-176]  and  thyroid  [177,  178],  is  resection  of  the  tumor  and  surrounding  tissue. 
Excised  tissue  is  evaluated  using  conventional  histological  examination  in  which  the  tissue  is 
fixed  in  formaldehyde,  dehydrated,  embedded  in  paraffin,  cut  into  micron  thin  slices, 
deparaffinized,  rehydrated,  stained,  dehydrated  and  finally  mounted  on  a  microscope  slide  for 
visual  inspection.  The  large  number  of  processing  steps  require  considerable  time,  usually  in 
excess  of  one  day  [179].  The  resulting  delay  between  surgery  and  pathological  assessment  poses 
complications  when  surgical  treatment  depends  on  the  results  of  histopathology.  For  example,  in 
the  case  of  early  breast  cancer,  breast-conserving  therapy  and  radiation  is  the  standard  of  care 
because  it  achieves  survival  and  local  recurrence  rates  equivalent  to  those  of  mastectomy  while 
providing  superior  cosmetic  outcomes,  reduced  morbidity  associated  with  lymphedema,  and 
improved  post-operative  quality  of  life  [180-184],  However,  local  recurrence  following  breast- 
conserving  surgery  is  strongly  correlated  with  margin  status  as  determined  by  post-operative 
histopathology  [185,  186],  Owing  to  the  limited  sensitivity  of  gross  examination,  residual  tumor 
is  present  near  the  surgical  margin  in  up  to  40%  of  cases  [187,  188].  In  the  event  of  a  positive 
margin  on  histopathology,  patients  typically  require  a  repeat  surgery  with  associated  morbidity 
and  delays  in  adjuvant  therapy.  Similarly,  in  the  treatment  of  papillary  thyroid  neoplasm,  a 


47 


finding  of  vascular  or  capsular  invasion  following  a  lobectomy  typically  requires  a  second 
surgery  to  perform  a  total  or  near-total  thyroidectomy  [178], 

Unfortunately,  the  long  processing  times  precludes  the  use  of  conventional  paraffin 
embedded  histopathology  for  the  intraoperative  evaluation  of  excised  tissue.  An  alternative, 
frozen  sections  analysis  (FSA),  replaces  dehydration  with  simple  freezing,  allowing  rapid 
intraoperative  evaluation,  but  has  limited  acceptance  owing  to  the  processing  time  per  section 
(-20  min)  and  low  sensitivity  and  specificity  compared  to  conventional  histology  [189-191]. 
Consequently,  relatively  few  surgeons  currently  use  FSA  for  intraoperative  evaluation  of  breast 
pathology  [192].  These  limitations  are  even  more  acute  in  the  case  of  papillary  thyroid 
neoplasms,  where  sensitivities  as  low  as  17%  using  FSA  have  been  reported  [193]  and  the  use  of 
FSA  for  papillary  thyroid  neoplasm  is  no  longer  recommended  [194].  There  is  a  need  for 
advanced  imaging  technologies  that  perfonn  rapid,  high  resolution  imaging  of  unprocessed  tissue 
samples  in  the  surgical  or  pathology  lab  setting. 

During  the  previous  reporting  period,  a  major  focus  of  our  research  was  the  development  of 
transfonnative  new  imaging  technologies  that  could  augment  or  eventually  replace  the  classic 
histopathology  techniques  in  applications  where  long  processing  times  interfere  with  patient 
treatment.  As  a  result  of  these  efforts,  we  developed  new  advanced  imaging  capabilities  based  on 
optical  coherence  microscopy  (OCM)  and  multiphoton  microscopy  (MPM)  that  enable  real-time 
visualization  of  cellular  structure  in  thick  sections  of  unprocessed  human  tissue.  These 
techniques  replace  the  lengthy  dehydration,  embedding,  microtoming  and  staining  process  with 
optical  sectioning,  which  enables  imaging  planes  corresponding  to  thin  histological  sections 
without  the  need  to  physically  cut  the  tissue.  These  methods  allow  rapid  evaluation  of  tissue 
pathology  in  a  clinical  setting  (e.g.  intrasurgical  imaging)  while  still  preserving  the  tissue  to 
enable  subsequent  conventional  histological  processing  (e.g.  post-operative  evaluation). 

6,1  Multiphoton  Microscopy  for  Rapid  Evaluation  of  Pathology 

Introduction  to  Multiphoton  Microscopy  (MPM) 

Multiphoton  microscopy  has  had  a  revolutionary  impact  on  fundamental  scientific  research 
and  biological  imaging  in  fields  ranging  from  cell  biology  to  neurosciences  [195,  196].  MPM  can 
achieve  cellular  resolution  and  high  contrast  axial  sectioning  in  thick  tissue  specimens  by  point 
scanning  short  pulse  near-infrared  light  across  tissue.  However,  MPM  can  also  provide 
molecular  specific  contrast  using  both  extrinsic  probes  and  intrinsic  fluorescence,  enabling  it  to 
visualize  many  of  the  same  biochemical  markers  that  are  the  hallmarks  of  pathology  in 
conventional  histological  examination.  The  excellent  resolution  and  molecular  contrast  available 
with  MPM  are  particularly  well-suited  for  visualizing  subcellular-resolution  features,  such  as 
nuclear  size  and  shape,  and  morphological  changes  in  tumor  extracellular  matrix  and  stroma, 
such  as  reorganization  of  collagen,  which  are  important  hallmarks  of  neoplasia  [197-200]. 

A  major  advantage  of  MPM  is  the  ability  to  use  fluorescent  probes  that  rapidly  target  specific 
cellular  agents,  including  DNA,  specific  proteins  and  cellular  receptors.  By  carefully  choosing 
contrast  agents,  MPM  images  can  be  generated  from  unfixed  tissue  that  closely  approximate 
more  familiar  microscopy  techniques,  such  as  the  use  of  DNA-specific  dyes  and  transmission 
wide-field  microscopy  in  conventional  clinical  histopathology,  but  without  the  lengthy  (20+  hour) 
processing.  This  concept  is  demonstrated  in  Figure  42,  where  an  intact,  unprocessed  human 
breast  tissue  specimen  has  been  stained  analogously  to  conventional  H&E  histology  and  then 
optically  sectioned  using  MPM.  Consequently,  MPM  has  the  potential  to  enable  real-time 
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evaluation  of  microscopic  pathology  in  a  surgical  setting  while  retaining  familiar  markers  of 
pathology  from  conventional  histology. 


Figure  42.  Virtual  3D  H&E  histopathology  performed  using  multiphoton  microscopy  (MPM)  of  freshly 
excised,  unfixed  normal  human  breast  tissue  showing  collagen  fibers  (red  through  pink)  and  cell  nuclei 
(purple). 

Preliminary  Multiphoton  Microscopy  Study  of  Breast  Surgical  Specimens 

Working  in  collaboration  with  Dr.  James  Connolly,  M.D.  at  the  Beth  Israel  Deaconess 
Medical  Center  and  Harvard  Medical  School,  we  investigated  the  sensitivity  of  MPM  by 
evaluating  freshly  excised,  human  breast  lumpectomy  specimens  for  the  presence  of  cancerous 
tissue  on  simulated  surgical  margins.  Discarded  (not  required  for  patient  diagnosis)  surgical 
specimens  were  rapidly  stained  with  a  nuclear  contrast  agent  (acridine  orange),  mounted  against 
a  coverslip,  and  imaged  using  two  spectral  channels,  corresponding  to  two-photon  fluorescence 
(TPF)  of  acridine  orange  to  visualize  nuclei  and  SHG  from  collagen  to  detect  stromal 
organization.  In  order  to  facilitate  clinical  interpretation  of  nonlinear  microscopy  images, 
pseudo-H&E  images  were  created  by  remapping  the  TPF  and  SHG  channels  to  a  colormap 
similar  to  that  of  standard  H&E  histology  (Figure  42).  This  use  of  two  channel  detection  and 
color  mapping  generates  tissue  contrast  that  is  analogous  to  standard  H&E  while  also  displaying 
images  using  analogous  colors.  This  method  has  the  powerful  advantage  that  the  MPM  images 
more  closely  match  the  interpretation  paradigm  used  by  pathologists  who  have  extensive  training 
on  H&E  images. 

In  order  to  enable  imaging  of  large  (several  square  centimeters)  tissue  areas  and  precise 
registration  with  conventional  H&E  histology,  large  field  of  view  images  were  generated  by 
rapidly  scanning  the  microscope  stage  and  specimen.  The  multiphoton  microscope  acquires  1024 
x  1024  pixel  images  of  a  -500  pm  x  500  pm  field  of  view  at  20  frames  per  second.  Multiple 
image  frames  were  mosaicked  in  post-processing  using  a  combination  of  unwarping  and 
nonlinear  intensity  blending  algorithms  (Figure  43).  The  resulting  composite  images  were  then 
converted  to  a  multilevel  tiled  pyramid  format  that  allowed  for  efficient  transmission, 
visualization,  zooming,  and  panning  of  large,  high  resolution  images  (Deep  Zoom,  Microsoft). 
This  automated  image  acquisition  and  processing  scheme  allowed  us  to  interrogate  entire  square 
centimeter  areas  of  breast  specimens  with  nuclear  resolution.  The  image  display  is  similar  to 
“Google  Maps”  and  enables  panning  a  large  field  of  view  to  rapidly  identify  regions  of  interest, 
then  zooming  to  a  high  magnification,  small  field  of  view  (Figure  44  (A)).  Images  at  multiple 
depth  levels  up  to  100  um  below  the  tissue  surface  were  also  obtained  to  generate  images 
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analogous  to  serial  sections  in  histology.  After  nonlinear  microscopy  imaging,  the  cassette 
containing  the  specimen  was  placed  in  10%  neutral  buffered  formalin  prior  to  processing, 
paraffin  embedding,  sectioning,  and  H&E  staining.  Corresponding  H&E  sections  were  scanned 
using  a  digital  slide  scanner  at  20x  magnification  (Aperio)  and  converted  to  the  same  file  format 
as  the  MPM  images  for  visualization  and  reading  (Figure  44  (B)). 


Figure  43.  Schematic  of  multiphoton  microscope  with  extended  field  of  view  and  color  remapping.  (A) 
Dual  detection  channels  are  used  to  collect  TPF  and  SHG  signals  for  nuclear  and  stromal  contrast, 
respectively.  A  1024  x  1024  pixel  image  can  be  acquired  at  a  rate  of  20  fps.  (B)  Overlapping  fields  are 
automatically  acquired  and  (C)  field  distortion  correction,  pseudo-H&E  color  remapping,  and  (D)  image 
mosaicking  is  applied  in  post-processing. 
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Figure  44.  Human  breast  specimen  with  ductal  carcinoma  in  situ  (DCIS).  (A)  MPM  image  and  (B) 
corresponding  H&E  histology  showing  different  subtypes  of  DCIS  (scale  bar:  1  mm).  Finger-like 
projections  of  abnormal  cells  invading  the  duct  lumen  are  indicative  of  micropapillary  DCIS  (red  box,  scale 
bar:  250  pm).  Holes  of  cellular  debris  and  comedo  plug  are  indicative  of  comedo  DCIS  (green  box,  scale 
bar:  250  pm).  Individual  cell  nuclei  at  comedo  DCIS  border  exhibit  prominent  nucleoli  consistent  with 
high-grade  DCIS  (blue  box,  scale  bar:  50  pm).  Field  of  view:  8.25  mm  x  5  mm  (660  MPixels). 

MPM  Data:  http://tinvurl.com/cwvly6s 
H&E  Data:  http://tinyurl.com/9teff9a 

We  assembled  a  dataset  of  168  specimens  from  50  patients  imaged  with  both  MPM  and 
conventional  H+E  histology.  The  dataset  was  divided  into  a  30  specimen  training  set  and  138 
reading  set  that  included  paired  MPM  and  H&E  images.  In  the  preliminary  study,  each  specimen 
was  diagnosed  based  on  the  following  tissue  subclassifications:  normal  breast  tissue, 
fibroadenoma,  fibrocystic  changes,  biopsy  site  changes,  flat  epithelial  atypia,  lobular  carcinoma 
in  situ  (LCIS),  DCIS,  invasive  ductal  carcinoma  (IDC),  and  invasive  lobular  carcinoma  (ILC). 
Similar  histological  subclassifications  were  then  grouped  during  data  analysis  to  determine  the 
diagnostic  accuracy  for  detecting  benign,  atypical,  and  malignant  tissues  using  nonlinear 
microscopy.  Three  pathologists  read  the  nonlinear  microscopy  and  H&E  images.  Each  reader 
was  initially  trained  on  paired  MPM  and  H&E  data  and  then  asked  to  read  the  138  nonlinear 
microscopy  images  while  blinded  to  the  results  of  H&E  histology.  The  three  pathologists  then 
read  the  H&E  images  one  week  later  to  reduce  memory  confounders.  The  blinded  reading  of 
MPM  images  achieved  an  average  95.4%  (C.I.  90.3  to  98.0)  sensitivity  and  93.3%  (C.I.  89.3  to 
95.9)  specificity,  respectively,  for  distinguishing  between  invasive  cancer  and  DCIS  versus 
benign  tissue  as  compared  with  the  corresponding  H&E  diagnosis.  The  negative  predictive  value 
was  97.1%  and  accuracy  was  94.1%. 

These  results  demonstrate  that  MPM  can  reproduce  the  diagnostic  capability  of  conventional 
histopathology  without  lengthy  fixation  or  processing.  These  results  were  published  in  in  PNAS, 
1 1 1,  43,  2014,  Tao,  et  al  [3].  The  results  from  this  study  enabled  us  to  develop  a  new  NIH  funded 
research  program  (R01  CA178636-02,  “Intraoperative  real  time  breast  cancer  margin  assessment 
with  nonlinear  microscopy”)  and  demonstrate  the  effectiveness  of  AFOSR  funding  research  in 
new  program  development.  In  this  new  program,  we  plan  to  develop  a  wide  field  of  view 
nonlinear  microscope  designed  for  the  pathology  laboratory  and  working  with  a  team  of 
pathologists,  radiologists  and  surgeons,  demonstrate  that  MPM  imaging  of  breast  surgical 
specimens  can  be  used  to  guide  surgical  resection  and  reduce  the  rate  of  second  surgeries. 

Development  of  Ultrahigh  Speed,  Wide  Field  of  View  Multiphoton  Microscopy  for  Pathology 

In  spite  of  its  high  sensitivity  and  specificity,  there  is  an  incompatibility  between  the  field  of 
view  of  MPM,  which  requires  high  numerical  aperture  imaging  to  generate  nonlinear  excitation, 
and  the  wide  field  of  view,  low  magnification  (1-4  mm)  typically  employed  to  evaluate  histology. 
In  our  pilot  study,  this  was  addressed  through  the  mosaicking  of  hundreds  of  rapidly  imaged  500 
pm  x  500  pm  fields  to  generate  a  composite  image  comparable  to  conventional,  low- 
magnification  histopathology.  However,  in  spite  of  the  relatively  high  imaging  speed,  this 
approach  alone  is  insufficient  owing  to  the  very  large  image  size  and  consequently  long  imaging 
times  required  to  generate  a  two  photon  dataset  comparable  to  a  conventional  pathology  slide. 

To  address  this  limitation,  we  are  collaborating  with  Thorlabs  Inc.’s  Advanced  Imaging 
Group  to  develop  next-generation  multiphoton  microscopy  optics  that  will  enable  wide  field  of 
view  over  areas  comparable  to  conventional  low  magnification  (5-1  Ox)  brightfield  microscopy 
but  preserving  high  numerical  aperture  beam  focusing  and  fluorescence  collection.  The 
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technology  we  are  developing  will  enable  rapid  (>2  fps)  imaging  over  a  4  mm  x  4  mm  field  of 
view  (16  mm  area)  in  unsectioned  human  surgical  specimens  using  a  combination  of  highly 
optimized  laser  scanning  optics  and  a  novel  objective  design  developed  by  Thorlabs.  With  this 
capability,  a  pathologist  will  be  able  to  directly  evaluate  surgical  specimens  without  the  need  for 
mosaicking.  Alternatively,  by  rapid  translation  of  the  sample,  we  estimate  that  the  system  will  be 
able  to  scan  and  mosaic  approximately  1800  mm2  per  minute,  enabling  rapid  intraoperative 
imaging  of  entire  surgical  specimens. 

Development  of  Gigavoxel  Multiphoton  Fluorescent  Lifetime  Imaging  Technology 

While  MPM  can  be  performed  with  extrinsic  contrast  agents  to  label  specific  cellular  or 
tissue  components  or  even  cellular  receptors,  these  agents  are  typically  not  approved  for  in  vivo 
applications.  For  applications  such  as  in  vivo  imaging,  where  extrinsic  contrast  is  not  pennissible, 
the  limited  intrinsic  contrast  of  tissue  may  not  be  sufficient  for  diagnostic  assessment. 
Fluorescent  lifetime  imaging  microscopy  (FLIM)  is  a  complementary  imaging  modality  to  MPM 
that  enhances  intrinsic  fluorescent  contrast.  FLIM  can  be  perfonned  using  MPM  for  axial 
sectioning  (MPM-FLIM)  [201].  In  FLIM,  single  or  two  photon  excitation  excites  a  fluorophore 
which  decays  with  a  characteristic  time  constant  that  depends  on  the  specific  fluorophore,  pH, 
temperature,  local  proteins  and  other  characteristics  [202-205].  FLIM  therefore  provides 
information  on  the  microenvironment  of  individual  cells  or  entire  tissues,  which  could  reflect 
metabolic  changes  associated  with  malignancy  [202].  MPM-FLIM  has  shown  promise  as  a 
means  of  detecting  pathology  in  human  breast  specimens  [206]  and  animal  breast  cancer  models 
[207]  by  imaging  changes  in  NADH  and  FAD  lifetimes  associated  with  tumor  metabolism  as 
well  as  changes  in  stroma. 

While  high  speed  imaging  and  mosaicking  of  samples  is  an  established  technique  for  wide- 
area  MPM  imaging  [208],  to  date  MPM-FLIM  has  been  limited  to  very  narrow  fields  of  view 
because  of  the  extremely  slow  acquisition  rates  of  existing  MPM-FLIM  systems  which  use  time 
correlated  single  photon  counting  (TCSPC)  [209-211],  It  should  be  noted  that  FLIM  without 
MPM  can  be  implemented  in  a  wide  field  geometry  significantly  faster  than  TCSPC  [212,  213], 
but  does  not  provide  the  depth  sectioning  required  for  imaging  intact  tissue  specimens  and  still 
suffers  from  relatively  low  imaging  speed  compared  to  MPM.  Therefore,  although  FLIM  has 
shown  promise  as  a  means  for  imaging  and  detecting  cancer  pathology,  it  has  not  been  possible 
perform  a  systematic  evaluation  of  MPM-FLIM  in  surgical  specimens  to  assess  its  potential  for 
intrasurgical  evaluation. 

The  limited  throughput  of  TCSPC  occurs  because  it  is  necessary  to  illuminate  individual 
pixels  at  sufficiently  low  power,  such  that  on  average  during  any  given  sampling  period,  only  a 
single  photon  arrives.  With  knowledge  of  the  instrument  response  function,  the  arrival  time  of  a 
single  photon  can  be  estimated  with  sub-nanosecond  accuracy.  However,  the  need  to  individually 
count  photons  greatly  limits  throughput  [209-211],  To  address  this  limitation,  we  have  developed 
a  fundamentally  new  approach  to  MPM-FLIM  measurement  based  on  optically  clocked,  direct 
temporal  sampling  of  fluorescent  emission.  In  this  method  (Figure  44),  adapted  from  optically 
clocked,  swept  source  OCT,  a  fast  digitizer  is  phase-locked  to  a  multiple  of  the  excitation  laser 
repetition  rate.  Individual  laser  pulses  are  imaged  onto  the  sample  at  nonnal  MPM  power  levels, 
generating  many  thousands  of  photons  per  excitation  and  allowing  measurement  of  the 
fluorescent  decay  time  at  rates  that  are  3-4  orders  of  magnitude  faster  than  TCSPC. 
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Figure  45.  Schematic  diagram  of  the  prototype  high  speed  MPM-FLIM  technology  showing  the  signal 
path  (upper  left  path)  and  the  reference  clock  (lower  right  clock).  A  GHz  sample  rate  A/D  combines  the 
sample  and  reference  signals  to  extract  time-resolved  fluorescence  without  single-photon  counting. 


To  explore  the  concept  of  optically  clocked  direct  temporal  sampling,  a  prototype  MPM- 
FLIM  system  was  developed  and  used  to  image  freshly  excised  human  breast  surgical  specimens 
(Figure  45).  Images  were  acquiring  with  a  single  spectral  channel  (375-700  ntn)  at  a  rate  of  10 
megapixels  per  second  over  a  1.1  x  0.8  mm  field  of  view  Figure  46  (A).  An  enlarged  view 
(Figure  46  (B))  shows  a  breast  ductal  structure  with  a  well-resolved  supporting  basement 
membrane  (yellow/green),  columnar  secretory  cells  (turquoise),  and  metabolically  inactive 
lumen  (violet).  This  preliminary  data  demonstrates  record  imaging  speeds  and  strongly  suggests 
that  TPM-FLIM  imaging  of  human  tissue  is  feasible  at  speeds  that  will  enable  surveying  of 
centimeter-scale  pathology  specimens. 


Figure  46.  16  megapixel  MPM-FLIM  image  of  freshly  excised  normal  human  breast  tissue  acaquired  at  10 
megapixels/s.  (A)  1.1  mm  wide  field  image  acquired  in  1.6  seconds.  (B)  Enlarged  view  of  a  ductal  unit 
showing  individual  cells  arranged  with  a  well-defined  lumen.  Scale  bars:  100  pm. 
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6.2  Development  of  Next  Generation  Optical  Coherence  Microscopy  Technology 

While  techniques  based  on  nonlinear  excitation  of  fluorescence  provide  molecular  contrast 
and  very  high  resolution,  the  need  for  tight  temporal  and  spatial  confinement  of  the  excitation 
beam  requires  high  NA  optics  and  are  strongly  affected  by  optical  aberration  and  scattering. 
This  limitation,  along  with  the  shorter,  more  strongly  attenuated  fluorescence  emission 
wavelength,  limit  the  maximum  imaging  depth  in  tissue.  Furthermore,  although  techniques  exist 
for  parallel  imaging  of  multiple  depths,  they  are  complex  and  impose  further  limitations  on 
maximum  imaging  depth.  Consequently,  MPM  is  most  effective  for  imaging  en  face  planes 
which  are  at  shallow  depths.  In  contrast,  interferometric  techniques  can  simultaneously  image 
many  depths  in  parallel,  operate  effectively  at  a  wide  range  of  numerical  apertures,  and  can 
image  at  longer  wavelengths  without  the  need  to  detect  short  (usually  visible)  wavelength 
fluorescence.  These  advantages  make  interferometric  techniques  effective  for  imaging  specimens 
that  have  three-dimensional  structure  or  specimens  for  which  the  region  of  interest  is  more  than 
100-200  um  below  the  surface. 

Optical  Coherence  Microscopy  (OCM) 

Optical  coherence  microscopy  (OCM)  combines  optical  coherence  with  confocal  microscopy 
to  enable  depth  resolved  visualization  of  specimens  with  cellular  resolution  [131,  214].  The 
addition  of  coherent  detection  greatly  increases  imaging  depth,  while  enabling  multiple  depths  to 
be  imaged  in  parallel,  enabling  higher  imaging  rates.  In  addition,  because  it  uses  interferometry, 
OCM  can  provide  depth-resolved,  quantitative  phase  information  from  materials  or  biological 
specimens,  which  offers  unique  advantages  over  other  phase  microscopy  techniques  such  as 
phase  contrast  and  differential  interferential  interference  contrast  microscopy  [215,  216]. 
Furthermore,  by  employing  polarization-sensitive  detection  techniques,  OCM  can  provide  depth 
resolved  birefringence  information  on  the  samples,  including  retardation,  axis  of  retardation  and 
degree  of  polarization  non-uniformity. 


Figure  47.  High  NA,  short  confocal  parameter  (a)  and  low  NA,  long  confocal  parameter  (b)  imaging 
regimes,  (c)  Dependence  of  the  transverse  spot  size  on  NA  of  the  objective,  f  :  coherence  length,  b: 
confocal  parameter. 

Figure  47  (A,  B)  shows  two  general  operating  regimes  for  OCM  imaging.  In  the  first  mode, 
a  relatively  low  NA  objective  is  employed,  which  results  in  a  long  confocal  parameter  for  the 
imaging  beam,  and  consequently  makes  it  possible  to  image  over  longer  depth  ranges.  However, 
an  increase  in  confocal  parameter  comes  at  a  cost  of  decreased  transverse  resolution.  In 
microscopic  imaging  applications  where  high  transverse  resolution  is  desired,  high  NA 
objectives  must  be  used.  As  seen  Figure  47  (B),  the  high  NA  configuration,  while  improving  the 
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transverse  resolution,  reduces  the  confocal  parameter  significantly.  OCM,  however,  circumvents 
this  trade-off  because  the  axial  resolution  is  not  set  by  the  NA  of  the  objective;  but  by  the 
coherence  length  or  axial  point  spread  function  of  the  coherent  detection.  Therefore,  the  NA  in 
OCM  is  selected  based  on  the  desired  spot  size.  As  shown  in  Figure  47  (C),  even  a  relatively 
low  NA  objective  gives  <3  um  transverse  resolution,  which  is  in  most  cases  sufficient  for 
achieving  cellular  resolution. 

Our  group  first  developed  and  demonstrated  OCM  in  1994  [131].  We  recently  demonstrated 
OCM  imaging  with  a  time  domain  system  at  a  speed  of  2  frames  per  second  (fps)  over  a  400  um 
x  400  um  FOV  [217].  This  system  was  subsequently  developed  into  an  integrated  3D-OCT  and 
OCM  instrument,  enabling  multiscale  imaging  of  cancer  pathology.  Working  in  collaboration 
with  Dr.  James  Connolly,  M.D.  at  the  pathology  department  of  Beth  Israel  Deaconess  Medical 
Center  (BIDMC)  and  Harvard  Medical  School,  we  performed  a  series  of  studies  using  integrated 
3D-OCT  and  OCM  to  image  nonnal  and  pathological  human  tissue  ex  vivo,  including 
gastrointestinal  tract,  thyroid,  breast,  and  kidney;  comparing  OCT/OCM  images  with  the  gold 
standard  histology  to  identify  the  characteristic  features  of  pathology  [135,  217-219].  During  the 
last  contract  period,  we  developed  advanced  technology  for  OCM  which  dramatically  improves 
image  acquisition  data  rates. 

High  Speed  Fourier  Domain  Optical  Coherence  Microscopy 

In  standard  microscopy,  features  are  viewed  in  the  en  face  plane  at  a  given  depth.  Therefore, 
although  Fourier  domain  OCT  methods  perform  superior  to  time  domain  due  to  the  ability  to 
image  many  depths  in  parallel,  if  the  image  is  to  be  acquired  from  only  a  single  depth,  this 
advantage  is  not  readily  appreciable.  In  addition,  since  each  pixel  in  an  en  face  image  requires  an 
A-scan,  extremely  high  imaging  speeds  are  required.  For  this  reason,  first  generation  Fourier 
domain  OCT  systems  were  not  suitable  for  OCM  and  this  delayed  the  development  of  Fourier 
domain  OCM  techniques.  However,  with  advances  in  light  source  and  detection  technology,  it  is 
now  possible  to  reach  axial  scan  rates  of  several  hundred  kHz  to  MHz.  Fourier  domain  OCM  has 
the  advantage  that  images  from  multiple  depths  are  acquired  simultaneously,  such  several  en  face 
planes  can  be  acquired  without  the  need  for  repeated  scanning  or  mechanically  translating  the 
interferometer  reference  arm.  In  addition  to  greatly  decreasing  the  time  required  to 
volumetrically  image  the  sample,  this  also  significantly  relaxes  the  need  to  match  the  coherence 
gate  depth  position  and  the  focal  plane  depth  position.  This  is  an  important  feature  for  integrating 
the  technology  with  confocal,  fluorescence  and  nonlinear  microscopes,  which  usually  have 
variations  of  path  length  with  scanning.  In  addition,  the  ability  to  measure  phase  of  the  back 
reflected  light  also  enables  numerical  algorithms  that  can  correct  for  aberrations  or  the 
defocusing  effects.  During  the  last  contract  period  we  investigated  the  development  of  Fourier 
domain  OCM  technology  in  both  swept  source  and  spectral  domain  configurations  for  pathology 
imaging. 

The  first  demonstration  of  swept  source  OCM  was  performed  by  our  group  in  2007  using  a 
high  speed  Fourier-domain  mode-locked  (FDML)  laser  [220].  OCM  imaging  was  performed  at 
42  kHz  axial  scan  rate,  with  imaging  resolutions  of  1.6  um  x  8  um  (transverse  x  axial)  over  a 
-220  um  x  220  um  field  of  view.  More  recently  our  collaboration  with  Praevium  /  Thorlabs  on 
the  development  of  high  speed  VCSEL  swept  source  technology  has  allowed  us  to  substantially 
improve  the  performance  over  earlier  swept  source  OCM  [221,  222].  The  VCSEL  laser  supports 
variable  sweep  rates  from  50  kHz  to  1.2  MHz  while  maintaining  a  sweep  range  of  >100  nm.  We 
have  demonstrated  swept  source  OCM  using  the  VCSEL  light  source  at  both  1050  nm  and  1310 
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nm  wavelengths,  with  an  axial  scan  rate  of  100  kHz  and  560  kHz,  respectively,  with  an 
adjustable  transverse  resolution  between  0.86  um  to  3.42  um  and  an  axial  resolution  of  8.1  um, 
comparable  to  conventional  histological  sections  [223]. 


Figure  48:  (A-E):  Depth  resolved  OCM  imaging  of  a  fresh  ex  vivo  human  thyroid  specimen,  with  a  10X 
objective,  (A)  Volume  rendering  emphasizing  that  arbitrary  planes  can  be  selected  for  visualization.  For  the 
imaging  planes  indicated  by  the  colored  lines,  reconstructed  cross-sectional  and  en  face  images  are  shown 
in  (B-E).  (C-E)  are  en  face  images  from  50  pm,  130  pm  and  180  upm  below  the  surface  of  the  specimen, 
respectively.  The  cross-sectional  image  in  (B)  is  displayed  using  logarithmic  scale,  whereas  (C-E)  are 
displayed  using  square  root  scale.  (G-J):  Ex-vivo  OCM  images  (G-I)  and  corresponding  histology  (J-L) 
from  fresh  normal  human  kidney  (G,  H,  J,  K)  and  clear  cell  renal  cell  carcinoma  specimens  (I,  L)  specimen 
obtained  using  40X  (G,  I)  and  20X  (H)  water  immersion  objectives.  (G,  I,  J,  L)  are  from  renal  cortex, 
whereas  (H,  K)  are  from  renal  medulla.  CT:  Convoluted  tubules,  G:  Glomerulus,  T:  Collecting  ducts. 
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Figure  48  shows  representative  images  taken  with  at  1310  nm  using  VCSEL  OCM  with  10X 
(Figure  48  (A-E)),  20X  (Figure  48  (G,  I))  and  40X  (Figure  48  (H))  objectives.  With  a  560  kHz 
axial  scan  rate  and  an  axial  scan  sampling  density  of  1024  x  1024  pixels,  an  imaging  rate  of  ~0.5 
volumes  per  second  was  achieved.  In  Figure  46  (A),  a  volumetric  rendering  of  an  OCM  dataset 
taken  from  freshly  excised  normal  human  thyroid  sample  is  shown.  From  this  dataset  one  can 
either  construct  images  in  the  en  face  direction,  similar  to  standard  microscopy  (Figure  48  (C- 
E)),  or  cross-sectional  images  with  very  high  transverse  resolution  (Figure  48  (B)).  Furthermore, 
it  is  possible  to  digitally  scan  in  depth  to  visualize  en  face  images  from  multiple  depths  (Figure 
48  (C-E)).  For  example,  Figure  48  (G-L)  shows  images  from  fresh  human  kidney  specimens, 
demonstrating  the  ability  of  SS-OCM  to  generate  structural  images  with  close  resemble  standard 
H&E  staining  in  fixed  specimens. 

Similar  results  can  also  be  achieved  using  spectral  domain  OCM.  Spectral  domain  detection 
has  the  advantage  that  extremely  broadband  light  sources  are  readily  available  around  800  nm, 
which  enables  very  high  axial  resolutions  (<5  pm)  that  is  not  yet  attainable  with  swept  source 
technology,  however  axial  scan  rates  are  limited  by  the  speed  of  line  scan  camera  technology. 
We  have  developed  a  SD-OCM  system  using  a  superluminescent  diode  (SLD)  light  source  that 
provides  a  4.2  um  axial  resolution  [224].  This  system  operated  at  a  210  kHz  axial  scan  rate, 
which  is  significantly  faster  than  most  spectral  domain  systems,  but  because  each  pixel  in  the  en 
face  image  requires  an  axial  scan,  the  en  face  frame  rates  were  limited. 


Figure  49.  Large  field  mosaic  OCM  image  of  a  nornial  human  kidney  specimen.  Image  is  constructed  by 
merging  30  frames  taken  with  the  40X/W  objective  producing  a  total  field  of  1.8  mm  x  2.1  mm.  (B)  shows 
a  zoomed  view  for  the  region  shown  with  dashed  lines  in  (A).  CT:  Convoluted  tubules,  G:  Glomerulus. 

Image  Tiling  and  Display  for  Large  Field  of  View 

OCM  at  high  numerical  aperture  has  a  limited  field  of  view  similar  to  MPM.  Even  using  a 
moderate  power  20X  objective  it  is  difficult  to  achieve  more  than  1  mm  FOV,  which  limits  the 
ability  to  evaluate  large  specimens  which  may  be  several  cm  in  dimension.  The  FOV  of  OCM 
can  be  increased  by  post-processing  and  tiling  multiple  smaller  fields  together  [223,  224]. 
Figure  49  shows  a  large  field  of  view,  high  resolution  OCM  image  of  human  kidney  constructed 
from  smaller  fields  that  have  been  tiled  together.  The  image  was  generated  by  combining  30 
frames  acquired  with  a  40X  objective  to  obtain  a  1.8  mm  by  2.1  mm  total  field.  Glomeruli  and 
convoluted  tubules  can  be  observed  throughout  the  imaging  field,  consistent  with  the 
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characteristics  of  normal  renal  cortex.  As  previously  described,  the  deep  zoom  display  fonnat 
allows  the  user  to  pan  around  the  image  and  zoom  to  regions  of  interest  enabling  a  more  detailed 
examination.  Figure  49  (B)  shows  an  example  where  a  region  with  a  glomerulus  surrounded  by 
convoluting  tubules  is  enlarged  to  show  finer  details. 

7.0  Other  Collaborations 

7.1  Computer-aided  Diagnosis  of  Plaque  Erosion 

The  two  most  frequent  underlying  mechanisms  for  sudden  cardiac  death  and  acute  coronary 
syndrome  (ACS)  are  plaque  rupture  and  plaque  erosion  [225-227].  Plaque  rupture  has  been  well 
characterized  both  ex  vivo  and  in  vivo  [225,  228],  However,  an  in  vivo  diagnosis  of  plaque 
erosion  in  patients  with  ACS  had  not  been  possible  due  to  the  lack  of  a  diagnostic  modality. 
Recently,  in  vivo  diagnosis  of  plaque  erosion  has  become  possible  using  intravascular  OCT  [6, 
229].  However,  the  diagnostic  algorithm  used  in  the  study  was  based  on  qualitative 
morphological  assessment,  which  has  potentially  high  inter-  and  intra-observer  variability. 
Plaque  erosion  has  distinct  pathologic  properties  including  proteoglycan-rich  and  smooth  muscle 
cell-rich  fibrous  regions  lacking  a  superficial  endothelial  layer  [226,  227,  230].  We  hypothesized 
that  quantitative  diagnostic  metrics  of  superficial  tissue  properties  can  provide  more  objective 
criteria,  which  may  enhance  in  vivo  diagnosis  of  plaque  erosion.  Furthermore,  quantitative 
metrics  can  allow  for  multiple  time  point  comparisons  of  the  same  lesion,  which  may  advance 
our  understanding  of  natural  history  of  plaque  erosion  and  its  response  to  treatment. 

Working  in  collaboration  with  Massachusetts  General  Hospital  (MGH)  OCT  Registry  led  by 
Dr.  Ik-Kyung  Jang,  we  have  developed  an  algorithm  for  in  vivo  quantitative  diagnosis  of  plaque 
erosion.  The  MGH  OCT  Registry  is  an  international  collaborative  effort  with  the  aim  to  facilitate 
clinical  intravascular  OCT  research  for  coronary  artery  diagnosis  and  treatment.  The  registry 
involves  more  than  20  international  sites,  with  a  clinical  intravascular  OCT  database  of  around 
3000  cases. 

A  total  of  42  patients  with  ACS  caused  by  plaque  erosion  were  included  in  the  study  [229]. 
Plaque  erosion  was  identified  according  to  the  previously  established  OCT  criteria.  Non-eroded 
fibrous  plaques  remote  from  the  erosion  site  within  the  same  vessel  were  used  as  controls 
(Figure  50).  Both  optical  properties  and  morphological  features  of  the  focal  eroded  region  as 
well  as  erosion-adjacent  region  were  quantified  and  compared  between  plaque  erosion  and 
control  fibrous  caps.  In  addition,  a  logistic  regression  model  was  developed  upon  selected 
features  for  automated  classification  of  plaque  erosion.  Compared  to  control  plaques,  eroded 
plaques  have  significantly  lower  surface  intensity  (p  <  0.00 1),  lower  region  of  interest  (ROI) 
intensity  (p  <  0.001),  lower  surface  normalized  standard  deviation  (NSD)  (p  <  0.001),  lower  ROI 
NSD  (p  <  0.001),  higher  optical  attenuation  (p  <  0.001),  larger  tissue  protrusion  area  (p  <  0.001) 
and  greater  surface  roughness  (p  <  0.001)  (Figure  51).  Erosion-adjacent  regions  also  have  lower 
ROI  NSD  (p  =  0.008),  higher  attenuation  (p  <  0.001),  and  greater  surface  roughness  (p  =  0.005). 
Using  the  logistic  regression  model  built  upon  the  quantitative  features,  plaque  erosion  can  be 
accurately  classified  against  intact  fibrous  plaques  (sensitivity:  100%,  specificity:  95.2%),  and 
the  likelihood  of  plaque  erosion  can  be  detennined  to  facilitate  diagnosis  (Figure  52).  The  entire 
vessel  can  also  be  automatically  labeled  in  a  continuous  colormap  indicating  the  likelihood  of 
plaque  erosion  (Figure  53). 
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Figure  50.  Representative  cases  of  plaque  erosion  and  a  control  fibrous  plaque.  Examples  of  plaque  erosion 
and  a  control  fibrous  plaque  included  for  analysis.  ROIs  were  traced  (blue  contour)  in  the  erosion-focal, 
erosion-adjacent  region,  and  in  a  control  fibrous  plaque  remote  (>  5  mm)  from  the  culprit  site. 
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Figure  51.  Quantitative  comparison  of  optical  and  morphological  properties  between  plaque  erosion-focal, 
erosion-adjacent,  and  control  fibrous  plaques.  Significant  optical  and  morphological  differences  are 
observed  in  all  the  metrics  between  erosion-focal  and  control  regions.  Erosion-adjacent  regions  also  have 
significantly  lower  ROI  NSD,  higher  attenuation  and  greater  surface  roughness  compared  with  control. 
Erosion- focal  regions  have  significantly  lower  surface  intensity,  lower  ROI  intensity,  lower  surface  NSD, 
lower  ROI  NSD,  larger  protrusion  area,  and  greater  surface  roughness  compared  with  erosion-adjacent 
region.  All  the  intensity  values  of  a  lesion  are  normalized  to  the  intensity  of  the  image  foreground. 


59 


Figure  52.  Classification  of  plaque  erosion  by  the  logistic  regression  model.  Examples  of  plaque  erosion 
(or  non-eroded  fibrous  plaques)  with  likelihood  ranging  from  100%  to  11%  as  determined  by  the  logistic 
regression  model.  The  likelihood  of  plaque  erosion  can  be  used  to  risk  stratify  the  lesion  progression. 
Definite  erosion  can  be  classified  by  simple  thresholding  (e.g.,  use  50%  as  the  cut-off  value). 
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Figure  53.  Automated  labeling  of  plaque  erosion.  Two  examples  of  OCT  images  are  labeled  automatically 
for  the  likelihood  of  plaque  erosion  using  the  proposed  algorithm.  At  each  angular  position,  the  region  of 
interest  (ROI)  was  selected  automatically  as  a  10  degree  region  bounded  radially  by  the  lumen  boundary 
and  the  vessel  wall.  The  algorithm  was  then  applied  to  the  ROI  and  the  likelihood  of  plaque  erosion  can  be 
generated  for  this  position.  The  entire  vessel  can  therefore  be  automatically  labeled  in  a  continuous 
colormap  indicating  the  likelihood  of  plaque  erosion. 

This  newly  developed  algorithm  promises  to  make  the  diagnosis  of  plaque  erosion  more 
accurate,  consistent,  and  simple.  With  the  ability  to  make  a  real  time  diagnosis,  cardiologists  will 
be  able  to  develop  treatment  strategies  tailored  for  each  individual  patient.  Future  studies  will  be 
focused  on  validating  the  algorithm  by  histopathology  studies,  and  applying  the  algorithm  to 
large  clinical  trials  to  correlate  the  prevalence  of  plaque  erosion  to  clinical  presentations. 
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7.2  OCT  for  Evaluation  of  Kidney  Transplantation 

There  is  currently  a  worldwide  shortage  of  kidneys  for  transplantation  due  mainly  to  the  fact 
that  there  is  no  reliable  means  to  determine  their  viability  (i.e.  kidneys  from  non-living  donors). 
Also,  following  transplantation,  kidneys  often  do  not  function  immediately  and  their  functional 
status  is  unclear.  In  studies  during  an  earlier  contract  period,  working  in  collaboration  with  Prof. 
Yu  Chen  at  University  of  Maryland  (formerly  a  postdoctoral  associate  in  our  group)  and  Prof. 
Peter  Andrews  at  Georgetown  University,  we  investigated  ultrahigh  resolution  OCT  imaging  of 
kidney  microstructure  in  the  rat  to  investigate  the  effect  of  renal  ischemia,  modeling  that  which 
occurs  in  kidney  transplantation.  3D-OCT  enabled  visualization  of  the  morphology  of  the 
uriniferous  tubules  and  the  renal  corpuscles.  OCT  revealed  dramatic  shrinkage  of  tubular  lumens 
due  to  swelling  of  the  lining  epithelium  during  ischemic  injury  and  3D  visualization  and 
volumetric  rendering  software  quantitatively  measured  volumetric  changes  in  tubular  lumens  in 
response  to  renal  ischemia.  Intravenous  infusion  of  mannitol  resulted  in  thinning  of  the  tubular 
walls  and  an  increase  in  the  tubular  lumen  diameters.  Mannitol  also  preventing  cell  swelling  that 
otherwise  resulted  in  shrinkage  of  proximal  tubule  lumens  during  ischemia.  These  studies 
showed  that  OCT  could  be  used  to  quantitatively  measure  structural  markers  which  were 
candidate  markers  for  kidney  viability  [231,  232]. 

Follow  on,  independent  studies  by  the  U.  Maryland  and  Georgetown  groups  investigated 
changes  in  transplant  protocol  for  living  donor  kidneys  to  reduce  ischemic  injury,  reducing  the 
time  of  mannitol  administration  from  30  to  15  min  or  less  before  clamping  the  renal  artery  [233]. 
OCT  imaging  showed  that  this  change  in  the  timing  of  mannitol  administration  protected  the 
human  donor  proximal  tubules  from  normothermic-induced  cell  swelling.  An  evaluation  of  post¬ 
transplant  recovery  of  renal  function  showed  that  patients  treated  with  this  modified  protocol 
returned  to  normal  renal  function  significantly  faster  than  those  treated  with  mannitol  30  min  or 
more  before  clamping  the  renal  artery.  The  U.  Maryland  and  Georgetown  groups  invited  us  to 
participate  in  an  NIH  program  (R01  DK094877-01,  “Non-invasive  evaluation  of  transplant 
kidney  using  OCT”,  P.I.  Peter  Andrews)  which  is  aimed  at  developing  and  applying  advanced 
OCT  imaging  technology  including  Doppler  OCT  and  OCTA  to  assess  structural  and  functional 
markers  of  kidney  function  in  transplant  surgery. 

We  and  our  collaborators  are  continuing  developing  advanced  ultrahigh  speed  swept  source 
OCT  technology  to  evaluate  the  structural  and  functional  markers  in  human  donor  kidneys  prior 
to  and  following  renal  transplant.  These  imaging  studies  will  be  correlated  with  post-transplant 
renal  function  in  order  to  develop  imaging  guidelines  and  associated  algorithms  for  OCT 
evaluation  of  kidney  viability  for  transplant.  The  proposed  studies  will  increase  the  number  of 
healthy  kidneys  available  for  transplantation  by  making  the  most  efficient  use  of  available  donor 
kidneys,  eliminating  the  possible  use  of  bad  donor  kidneys  (i.e.,  excessively  damaged  ischemic 
or  diseased  kidneys),  providing  a  measure  of  expected  post-transplant  renal  function,  and 
allowing  better  distinction  between  post-transplant  immunological  rejection  and  ischemic- 
induced  acute  renal  failure. 

Although  the  initial  research  leading  to  this  program  was  performed  in  2007,  this  is  an 
example  of  the  longer  term  impact  of  AFOSR  funding  on  clinical  care  as  well  as  new  program 
development. 
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